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SUMMARY 


This report documents the real-time algorithm used to 
generate the necessary video interface signals, for each of the 
selected oil rig/helicopter relative geometry scenarios, to drive 
an airborne color radar display. Typical display results are 
documented. Unique features of the real-time algorithm include 
(1) a circular list approach to target sorting and ray 
generation, (2) target shape and multiple target merging effect 
generation, and (3) sea clutter generation. 

Gratitude is expressed to Mr. George Clary of the NASA-Ames 
Research Center for his efforts on the project. In particular, 
his technical guidance on the radar math model and his work to 
provide in-flight radar video data for simulation validation was 
much appreciated. 
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I. INTRODUCTION 


1,1 BACKGROUND 

This report documents software developed by Systems Control 
Technology, Inc. (SCT) to provide a real-time simulation of an 
airborne radar for overwater approaches to oil rig platforms. 

This software was developed for NASA-Ames Research Center as part 
of the Rotorcraft All-Weather Operations Research Program. The 
simulation is being used to study advanced concepts for 
enhancement of airborne radar approaches (ARA) in order to reduce 
crew workload, improve approach tracking precision, and reduce 
weather minimums. ARA's are currently used for off-shore 
helicopter operations to and from oil rigs. 

With this purpose in mind, NASA-Ames needed a weather radar 
simulator for use in an existing fixed-base helicopter 
simulator. They awarded a contract to Computer Avionics 
Corporation (CAC) to develop the necessary hardware interface 
between the Sperry/RCA Primus 500 color weather radar and a Xerox 
Data Systems Sigma 9 simulation computer currently used by 
NASA-Ames. The radar interface hardware provides for two-way 
communication between the Sigma 9 and the radar display. 
Information obtained from the radar display through the interface 
includes antenna position, sweep direction, range setting, mode 
settings (weather/map and radar/beacon/both), radar and beacon 
gain settings, tilt, and sector scan (120 or 60 degrees). The 
radar interface receives from the simulation software encoded 
data defining the current target (oil rig) to be displayed. The 
ray definition consists of a sequence of intensity levels and 
ranges. The range defines where the radar is to display the new 
intensity level. 
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Figure 1.1 Overview of Simulation 
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1.2 OBJECTIVES 



SCT was awarded a contract with the objective of developing 
the software necessary to simulate overwater radar displays. The 
real-time simulation requirement for the software was a 
constraint requiring trade-off studies to achieve a reasonable 
degree of simulation fidelity using simplified math models and 
computationally efficient algorithms. Also, due to the projected 
research uses for the software modular programming, the use of a 
high-level language (FORTRAN) was required to achieve research 
flexibility. 

1.3 APPROACH 

An overview block diagram of the radar simulation software 
is depicted in Figure 1.1. The software consists of three basic 
sections. The first section is the initialization section which 
allows the operator to select a block of targets and to execute 
this code to define all necessary parameters before actual 
simulation. This code is only executed once before the 
simulation starts. The second block of code deals primarily with 
target manipulation. It is responsible for maintaining the 
target list in a circular linked list, maintaining pointers into 
the target list for targets in view of the radar and deciding 
what targets should be updated. The third area of code is 
responsible for modeling visible radar targets and sea clutter by 
inserting intensity level changes into a linked list. At the end 
of each pass, this list is then used to merge all targets and 
output the sequence of intensity levels and ranges to the radar 
interface . 

There are six subroutines in the radar simulation software. 
The first two deal only with initialization. 
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INIRDR 


Called by the simulation executive to initialize the 
radar simulation. Contains the X, Y locations of all 
targets to be simulated. Operator must define NBRTGT 
to indicate the first target in the block and NUMTGT to 
indicate the number of targets before INIRDR is called. 

INSTGT - Called by INIRDR to install a particular target into 
the linked list of targets. Used only for 
initialization. 


The remaining four subroutines are used when the simulation 
is running, and they perform the target manipulation and display 
generation tasks. 


CREATE - 


INSERT - 


UPDATE - 


HDGDST - 


The main subroutine called by the simulation executive 
to generate a radar vector. It inputs data from the 
radar interface and determines the antenna heading. It 
then processes targets starting with a target known to 
be behind the antenna sweep and continuing until it 
encounters a target beyond the antenna sweep. For each 
target visible, it generates a set of intensity level 
changes and inserts them into a linked list by a call 
to subroutine INSERT. It also determines which target 
should be updated with new relative position 
calculations and updates it with a call to subroutine 
UPDATE. It also inserts any intensity level changes 
caused by the sea clutter model. Finally, it uses the 
linked list of intensity level changes to create the 
encoded data for the radar display. 

Insert intensity level changes into the linked list. 
Each call inserts an increased level change at the two 
ranges supplied with the call. 

Updates a radar target position by a call to HDGDST. 

It then makes sure that the target is still in its 
correct position in the circular linked list. If not, 
it is moved to its correct position. 

Calculates the relative bearing, relative tilt and 
distance for a target. It is called by UPDATE and also 
INSTGT to calculate the target position. 


The constraint for real-time 
simulation software to execute in 
each radar ray generated. It was 
cycle time under 45 milliseconds; 
visual system and instrument panel 


operation requires the radar 
13 milliseconds or less for 
desired to keep the simulation 
the helicopter simulation with 
outputs requires 32 
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milliseconds per cycle, leaving up to 13 milliseconds of cycle 
time for the radar software. Using the 45 millisecond cycle 
time, the radr antenna moves 1-1/4 to 1-1/2 degrees per cycle. 
Therefore, each sof tware-generted ray is repeated until a new one 
is sent, and the software resolution is 1-1/4 to 1-1/2 degrees. 
This has not limited the display realism since the weather radar 
resolution is 6 to 8 degrees beamwidth. The information used to 
model the radar returns from point targets, beacons, and the sea 
was obtained from George Clary of NASA-Ames Research Center. 

Also, information regarding modeling of the threshold levels for 
the Sperry/RCA Primus 500 color radar and validation of the 
actual display were obtained with his assistance. Information 
regarding the radar theory and general description of the RCA 
Primus 500 radar is being documented in a NASA Technical Memo, 
"Simulation of a Weather Radar Display for Overwater Airborne 
Radar Approaches." 

The rest of this report is organized as follows. Chapter II 
reviews the key problems and issues of the radar simulation and 
describes the selected approach. Chapter III documents the 
high-level flow of the radar simulation. Chapter IV shows the 
radar display output from the simulation for various radar 
parameters and attempts to show typical results from varying one 
parameter at a time. Chapter V presents concluding remarks. 
Appendix A contains a complete source listing of the simulation 
software as well as documentation for each subroutine describing 
key variables and common areas. 



5 



This Page Intentionally Left Blank 



II. PROBLEM DISCUSSION AND APPROACH 


This section discusses the key problems and issues and then 
proceeds to describe the specific approach used in resolving each 
of them within the scope of this effort. 

2.1 OVERVIEW OF KEY PROBLEMS AND ISSUES 

Several major problems had to be resolved in order to meet 
the requirements of the real-time radar simulation^ These 
included the following. 

(1) Devising a method of maintaining the lists of targets 
in sorted order so as to minimize the number of targets processed 
per pass. This was deemed necessary because of the real-time 
requirements . 

(2) Determining a method to update a target's position in a 
timely manner in order to give an accurate position and to 
prevent the target from being updated while it is being 
displayed. This will prevent a target from being split. 

(3) Determining a time-efficient method for modeling the 
target shape. This method should accurately account for known 
radar theory. It must also reflect changes in the target range, 
range setting, radar mode, gain settings, beacon targets, and 
antenna tilt. 

(4) Devising a method to model sea clutter. Since little 
is known about the effects of sea clutter, an experimental model 
should be developed which will be time-efficient and will model 
some of the known or anticipated effects of sea clutter. It 
should take into account both the velocity and direction of the 
wind, thereby defining the "sea state." 

(5) Implementing a simple method to merge multiple targets 
and/or sea clutter. Since each target will be independently 
processed and the output of this processing will consist of 
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intensity changes at calculated ranges, then a sorted list of 
these intensity changes will need to be maintained. Before 
output to the actual hardware, those targets which overlap in 
range will have to be merged. 

2.2 TARGET TRACKING 

To satisfy requirements (1) and (2), it was decided to 
maintain the targets in a circular linked list in order of 
magnetic bearing from the helicopter's current position. Since 
the antenna sweeps in both the clockwise and counterclockwise 
direction, two pointers for each target were maintained. One is 
used for the clockwise direction and the other for the 
counterclockwise direction. Magnetic bearing was used as opposed 
to heading relative to the helicopter since magnetic bearing 
changes slowly with the helicopter's position while relative 
heading can change quickly as the helicopter changes heading. 
These pointers are stored in the TGTPTR array. 

In order to process only the targets pointed to by the 
antenna, two angular areas were defined. Both areas are defined 
relative to the current antenna direction. The first is the 
target area which defines the area in which targets may be' 
visible and should be processed. It is defined by the TGTWDT 
parameter which specifies the maximum difference in direction 
allowed between the target and the antenna direction. It is 
currently about 12 degrees, and therefore the target area extends 
from the antenna heading minus 12 degrees to the antenna heading 
plus 12 degrees. Once a target is in the target area, it can no 
longer have its position updated. This prevents the target from 
splitting since its position remains constant throughout the 
sweep past the target direction. 

The second angular area is the look-ahead area. It is 
defined in order to assure that the target's positions are up to 
date before entering the target area. This area extends from the 



forward edge of the target area about 6 degrees. This angle is 
defined by the LOKAHD parameter, but the actual size of the 
look-ahead area is determined by LOKAHD - TGTWDT. This area is 
always in front of the antenna sweep direction. When the antenna 
changes direction, the look-ahead area changes to the other side 
of the target area. When a target enters the look-ahead area, it 
becomes a primary candidate for updating its position. The 
proper size of the look-ahead area is determined by the number of 
targets that may enter the look-ahead area at the same time. 

Once a target enters the look-ahead area, it has about five 
passes to update its position before it enters the target area 
and becomes "frozen." If six targets came into the look-ahead 
area at once, then only five would have a chance to be updated 
before they entered the target area. The current value of LOKAHD 
seems appropriate since targets never enter the look-ahead area 
at that rate and the target's position changes very little in .6 
seconds; this is the time from when the target gets updated to 
the time the antenna sweeps by the target. A diagram of the 
target and look-ahead areas is contained in Figure 2.1. 

Targets have their position updated by two different 
priorities. Only one target is updated on each pass of the 
algorithm. First priority goes to targets that have entered into 
the look-ahead area. If no targets have entered the look-ahead 
area, a target is picked that is outside both the target and 
look-ahead areas. These secondary priority targets are selected 
by traversing clockwise through the target list. When the target 
area is encountered, it is skipped until a target is found past 
the target area. The actual updating is done by subroutine 
UPDATE. This routine updates the true heading and distance from 
the helicopter to the target. The distance computed is actually 
the slant range computed from the differences in X position, Y 
position, and altitude between the helicopter and the target. It 
also updates the tilt angle to the target and the log 10 value of 
the target's distance. If the target heading has changed so that 





TGTWDT = .2 RADIANS = 11.46 DEGREES 
LOKAHD = .3 RADIANS = 17.19 DEGREES 


Figure 2.1 Target and Look-Ahead Areas 
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it no longer lies between its clockwise and counterclockwise 
neighbor, the target pointers are changed by removing the updated 
target, finding its new position in the linked list, and then 
inserting it. 

2.3 TARGET MODELING AND SIMULATION 

Targets are simulated only if they are in the current target 
area. Then target modeling starts by computing the signal 
strength returned from the current target. The equations for 
calculating signal strength are derived from basic radar theory. 
The NASA technical memo, which is currently in work, will contain 
a detailed discussion of the radar theory involved in the 
simulation as well as discussions on the radar-specific values 
used in the radar equation. This memo will define the returned 
signal strength equatons for both normal and beacon targets, and 
the threshold levels for all three intensity levels. 

The basic radar equation for signal strength returned from a 
non-beacon oil rig target is: 

Return Power (dB) = Transmit Power - 30 * log (4*Pi) 

- 40 * log (Range) + 2 * Antenna Gain 

(in dB) 

+ 20 * log (Lambda) + 10 * log (Sigma) 

- TR Loss 

Several of these values are constants for purposes of this 
simulation. These include 

Transmit Power = 7000 watts = 68.45 dbm 

Lambda = .0319 meters = radar wavelength 

Sigma = 1000 square meters = target reflectivity 

TR Loss = 3.6 dB = SH-3 transmit/receive losses 

The antenna gain is computed from the maximum antenna gain minus 
antenna pointing error loss. The maximum antenna gain depends on 
the antenna size and can be computed by 6.4 + 200 log (antenna 
diameter in inches). For an 18" antenna, this gives a value of 
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31.5 db. The pointing error is calculated from the direction 
error and the tilt error. The gain loss is then calculated using 
the angular antenna error as an index into a table of antenna 
losses. The range is obtained from the target position 
information which was computed during the last update of the 
target. The range value is in nautical miles while the radar 
equation expects it in meters. The log of the conversion factor 
is therefore added in as a constant. The final radar equation 
for return power after all constants have been removed is 

Return Power (dB) = - 95.13 - 40* log (Range in NM) 

+ 2.00* antenna gain - TR loss 

The returned signal strength is then compared with the 
threshold levels. These define the minimum level for which the 
target is visible. The actual threshold levels depend on several 
factors which include pulse width, radar mode (map or weather), 
target distance and maximum antenna gain. Figures 2.2 and 2.3 
show the threshold levels for weather and map mode, respectively, 
as obtained from Sperry/RCA avionics for the RCA Primus 500 
radar. The MDS (minimum detectable signal) depends on the pulse 
width. For weather mode or map mode with a range of 50 miles or 
greater a long pulse width (2.35 s) is used and the MDS is 
lower. This is also used if the radar is in beacon or both 
modes. For map mode less than 50 miles in range, the shorter 
pulse width (0.6 s) is used; this gives a higher MDS. 

The 28 dB maximum antenna gain is the correct value for the 
12" antenna. If an 18" ^antenna is installed, the maximum antenna 
gain increases to 31.5 dB. This causes the returned signal 
strength to be increased by 7 dB. To correct for this, the 
threshold levels for the simulated radar are also increased 
whereas an actual hardware adjustment to the radar receiver is 
required on an aircraft. 

The target is modeled independently for each color 
intensity. The threshold level for the color intensity being 
modeled is subtracted from the returned signal strength. If this 
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yields a positive dB level, then the target is visible at that 
intensity. The length of the target for this intensity level is 
obtained from parameters defining the minimum length and a factor 
multipled by the dB level. The actual length is a function of 
the range selected, the color intensity, and the pulse width. 
Parameters also define the percentage of the target ahead and 
behind the actual target distance. Target shape and the 
parameters that define it are given in Figure 2.4. 

Targets with beacons on them are also modeled. Return power 
from beacons is defined as 

Return power (dB) = beacon transmit power 

- 20 log (4*PI) 

- 20 log (range in meters) 

+ radar antenna gain 

+ beacon antenna gain 

+ 20 log (lambda) 

- one-way loss 

Several of these values are constants for purposes of this 
simulation. These include: 

Beacon transmit power = 400 watts = 56.02 dBm 

Radar beacon antenna gain = 5 dB 
’Lambda = 0.0322 meters 

The range was again converted to nautical miles and the radar 
antenna gain is computed as before. This yields the return power 
as 

Return power (dB) = -56.16 - 20 * log (range in nm) 

+ antenna gain - one-way loss 

The threshold level for beacon returns is modeled as in 
Figure 2.5. There is no MDS defined. The dB level is obtained 
by subtracting the threshold level from the return power. If it 
is positive then the beacon target is visible and is modeled as a 
fixed length target. This length is the same as TGTLEN described 
in Figure 2.4. The beacon target appears at the maximum 
intensity while other non-beacon targets are prevented from 
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(MINPCT + PCTDB * OB) 
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ACTUAL PARAMETERS 


RANGE SETTING 

TARGET LENGTH 

PULSE WIDTH 

FACTOR 

2.5 NH 

.15 NM 

2.35 US 

2.0 

5.0 NM 

.15 NM 

0.6 iis 

1 .0 

10.0 NM 

.30 NM 



25.0 NM 

.40 NM 



50.0 NM 

.80 NM 

INTENSITY 

FACTOR 

100.0 NM 

1 .60 NM 

GREEN (LOW) 

1.0 



YELLOW (MED) 

0.65 



RED (HIGH) 

0.25 


MINPCT = MINIMUM LENGTH (DB=0) = 0.7 

PCTDB = ADDITIONAL LENGTH PER DB = 0.015 

ASTPCT = PERCENT OF TARGET BEHIND ACTUAL DISTANCE = 0.85 

AHDPCT = PERCENT OF TARGET AHEAD OF ACTUAL DISTANCE = 0.15 


Figure 2.4 Target Modeling Parameters 
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reaching maximum intensity. Because of the time delay between 
the beacon receiving an interrogation signal and the beam 
responding, the beacon target appears 0.6 miles behind the actual 
target. Thus, when the display mode is "both" (radar and beacon 
targets), the radar target with the beacon will appear at the 
correct distance, while the beacon target will appear 0.6 miles 
behind it. 

2.4 INITIALIZATION 

Several requirements were necessary for the initialization 
task. Initialization occurs only once prior to executing the 
simulation when subroutine INIRDR is executed. Many different 
target patterns were stored in the 100 loctions reserved for 
target descriptions. The operator must set two parameters (first 
target to be displayed and the total number of targets). The 
initialization procedure then calculates the position and bearing 
and then inserts it into the target pointers which maintain the 
clockwise and counterclockwise lists. This is repeated for the 
total number of targets. This procedure provides a means for the 
operator to define new target locations interactively and then to 
display them without the need for recompilation. 

2.5 SEA CLUTTER MODEL 

The sea clutter model is intended to model the effects of 
radar return from the sea. The sea clutter effect on actual 
radar is very complicated. There can exist very random returns 
from the sea with the amount of clutter ranging from almost solid 
with a few holes to an almost clear screen with a few solid 
patches. Because of the execution time constraints, several 
simplifying constraints were made. Once such constraint is that 
no holes would exist in the sea clutter. This allowed sea 
clutter to be defined as a single length of sea "noise." 
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Furthermore, it was decided that the sea clutter would only exist 
in the lowest intensity level. A separate pass through the 
algorithm would have to be made to generate a different intensity 
level and it was decided that the amount of actual sea clutter at 
a higher intensity level was not sufficient to justify more 
execution time. 

Therefore, the sea clutter model consists of calculating a 
return level from the sea and comparing it with the threshold 
level. The threshold level is calculated as in the radar target 
with the range being defined as the slant range between the radar 
antenna and the sea. The level of return from the sea accounts 
for several factors, which include the following. 

(1) The increasing reflectance of the sea as tilt 
decreases (pointing more directly at the sea); 

('2) The increasing reflectance as sea state 

increases. The sea state was defined in terms of 
wind velocity for purposes of this simulation. 

(3) The changes in reflectance a a function of heading 
relative to the sea. 

Other factors affecting the sea return include radar gain, pulse 
width, slant range, antenna gain, and miscellaneous transmission 
losses. 

After the sea return level is calculated, it is then 
compared with the threshold level. If the return level is 
greater than the threshold level, then an area of sea clutter is 
generated. This area is centered around the sea slant range. To 
provide more realism, a random noise factor was inserted in the 
calculation of sea clutter starting and ending ranges. 

Further discussion on the derivation of the sea clutter 
model is to be included in a NASA Technical Memo to be published 
later this year. 
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2.6 TARGET MERGING 


The final problem remaining was how to merge the 
independently generated targets when two or more targets 
overlap. The solution was to maintain a range sorted list of 
intensity changes. The intensity change defines both the 
previous intensity level and the new intensity level traversed. 
At every range, a count of each intensity level is maintained. 
If the intensity level produced is the same as the current 
intensity, then no change is necessary and the point is 
eliminated. An example of the merging of two targets is in 
Figure 2.6. 


The rules for merging intensity are as follows: 


2 level 1 

3-5 level 1 

6 or more level 1 

1 level 2 + 3 or more level 1 

1 level 2 + 1 or more level 2 


= level 1 (low intensity) 

= level 2 (medium intensity) 
= level 3 (high intensity) 

= level 3 
= level 3 


These rules were designed to 
large ranges (multiple targets in 
merging of targets at short range 


produce realistic merging 
close proximity) and also 
and high intensities. 


at 
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Figure 2.6 Target Merging Effects 
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III. PROGRAM DOCUMENTATION 


The methodology used in developing the real-time software 
was first to implement it on a system consisting of a PDF 11/34 
computer driving a Megatek graphics display which, in turn, was 
used to simulate the RCA radar. This allowed substantial 
development and debugging to be performed without access to the 
heavily used NASA simulation computer. After a baseline version 
had been developed and optimized for real-time execution, the 
program was converted to the NASA Sigma 9 computer. After 
several revisions of this baseline code on both the above- 
mentioned computers the final software, documented below, 
resulted. The work on the Sigma 9 facility was supported by 
Computer Sciences Corporation. Consequently, on-site personnel 
are acquainted with the code developed at a working level and are 
qualified to make future revisions. 

The following three Figures 3.1, 3.2 and 3.3 show the high 
level flow of the CREATE subroutine including target processing, 
modeling and updating. Further documentation including program 
description, listings and a list of variables is contained in 
Appendix A. This is the current software existing at NASA-Ames 
Research Center on the Sigma 9. 




Figure 3.1 Target Processing 
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Figure 3.2 Target Modeling 
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IV. DISPLAY OUTPUT RESULTS 


The objective of this effort was to produce visually 
realistic displays for typical overwater approaches to oil rig 
clusters and platforms. In order to allow project engineers to 
perform this activity it was necessary to (1) review currently 
available data on such approaches and (2) characterize the 
quality of the displays by providing hard copy examples of the 
display outputs for typical sets of key system parameters. 

As the work progressed, it turned out that very few flight 
test data sources were available; in fact, NASA-Ames project 
engineers appeared to have the best available visual description 
of the target returns. A flight was made in one of the NASA 
aircraft specifically for obtaining visual display data from oil 
rigs, and the NASA TM in preparation will show comparisons 
between photos of the simulation display and airborne display. 

As a mechanism to characterize the quality of the display 
generated, a series of runs wre made at the SCT facility and the 
results were drawn on a plotter. Table 4.1 summarizes the 
various runs. It should be noted that these runs were generated 
with SCT's software version of the radar simulation. The Sigma 9 
listing in Appendix A has modified some of the parameters to 
improve target shaping and sea clutter models. Therefore, these 
plots are not exactly as in the Sigma 9 version although the 
basic characteristics of the simulation are accurately 
represented . 

Figure 4.1 details the parameter values used for the 
simulation runs shown in Figures 4.2 through 4.45. Several 
observations should be made about the trial runs. Figures 4.2 to 
4.7 show the effects of target merging as the target distance 
goes from 1.75 to 70 miles. Also apparent is the pulse width 
change as the range in map mode is changed from 25 to 50 miles. 

As the pulse width changes from .6 us to 2.35 us the target 
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Table 4.1 

Display Characterization Cases 


FIGURE NO. 

COMMENTS 



4.2 - 4.7 

Effects of range on MAP mode 



4.8 - 4.12 

Effects of range on WEATHER mode 



4.13 

Both radar and beacon 



4.14 

Beacon only 



4.15 - 4.18 

Effects of radar gain 



4.19 - 4.22 

Effect of beacon gain 



4.23 - 4.29 

Effects of tilt 0-knot wind 

2000 ft 

4.30 - 4.34 

Effects of tilt 10-knot wind 

2000 

ft 

4.35 - 4.39 

Effects of tilt 10-knot wind 

2000 

ft 

4.40 - 4.45 

Effects of tilt 10- 30-knot wind 

500 

ft 
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FIGURE 
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0 
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length increases. Figures 4.8 to 4.12’ show the same runs in 
weather mode. The major difference between map and weather mode 
is the use of the 2.35 ps pulse width at all ranges in the 
weather mode. This causes the targets in weather mode to be 
thicker at ranges less than 50 miles than targets in map mode. 
Figure 4.13 shows the effect of placing the radar in both mode 
which displays both beacon and radar returns. All radar returns 
are suppressed from the highest intensity (red) while the beacon 
return is a constant thickness at the highest intensity. The 
beacon delay of .6 nautical miles places the beacon return behind 
the primary target. Figure 4.14 shows the beacon only display. 
Figures 4.15 to 4.18 shows the effects of the radar gain. As the 
gain voltage is lowered from 5 to 2 volts the targets shrink. 

This represents a change in gain from +3 to -21 dB compared to 
the preset gain equal to 4.625 volts. The targets disappear 
completely at radar gains of less than 2 volts. 

Figures 4.19 to 4.22 show similar effects on the beacon 
target when the beacon gain is reduced from 5 volts to 2 volts. 
Figures 4.23 to 4.29 show the effects of tilt on target 
appearance. When the radar tilt deviates from the optimal the 
target size diminishes. Also the first effects of sea clutter is 
shown in Figures 4.28 and 4.29. As the tilt decreases the 
angular width of the sea clutter increases while the range of the 
clutter decreases. Figures 4.30 to 4.39 show sea clutter effects 
as the wind velocity increases to 10 and 20 knots with tilts 
ranging from -1 to -12 degrees. Maximum range of the sea clutter 
is with a wind velocity of 20 knots and -1 degree tilt. The sea 
clutter in this case extends up to 50 miles. The last set of 
figures are sea clutter effects at 500 feet altitude. This gives 
shorter ranges of the sea clutter. Figure 4.45 shows sea clutter 
when the wind direction is perpendicular to the flight path. 

This causes the waves to be parallel to the flight path and less 
sea clutter return. The width of the sea clutter is maximum at 
the edges and minimal straight ahead as opposed to the other way 
around with a wind direction of 360 degrees. 
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A few typical cases of interest were also documented by 
photographing the RCA color display as simulated at NASA Ames 
Research Center. These photos are shown in Figures 4.46 through 
4.53. 
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Figure 4.46 Map Mode for Three-Target Situation 



Figure 4.47 Weather Mode for Three-Target Situation 
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Figure 4.48 Map Mode of Multiple Targets 



Figure 4.49 Weather Mode of Multiple Targets 
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Figure 4.50 Two-and-a-Half Mile Range Under 
Wind Conditions 



Ten Mile Range Under Wind Conditions 
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Figure 4.51 
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V. CONCLUDING REMARKS 


This report documents the software developed to drive a 
color radar display for overwater helicopter approaches to 
off-shore oil rig platforms. A variety of off-shore scenarios 
representing airborne radar approach conditions were defined in a 
software library and snapshots of the resulting output for 
interesting situations were documented in the report. 


Unique features of the real 
circular list approach to target 
target shape and multiple target 
(3) sea clutter generation. 


-time algorithm are (1) a 
sorting and ray generation, (2 
merging effect generation and 


) 
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APPENDIX A 


HELICOPTER RADAR DISPLAY SIMULATION PROGRAM LISTING 
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1 


SUBROUT IIMt CREATE 


DE3C R I P T I OhJ : 

THIS PROGRAM WAS WRITTEN FOR TmASA-AMES TO SIMULATE 
RADAR RETURNS IN A RCA PRIMUS 500 RADAR SCREEN. 

THIS SUBROUTINE CREATES A RADAR VECTOR IN A DIGITAL 
ENCODED FORM. THIS VECTOR IS IN A FORM COMPATIBLE 
WITH THE RADAR TARGET SIMULATOR WHICH DRIVES A PRIMUS 
500 RADAR SCOPE. THIS PROGRAM CREATES THESE VECTORS 
FOR SEA MODE ONLY. 

CRE.ATE IS CALLED EACH TIME A NEW VECTOR IS TO BE 
PRODUCED. IT INPUTS DATA FROM THE RADAR HARDWARE TO 
DETERMINE THE ANTENNA POSITION. DIRECTION OF TRAVEL 
AND RANGE SETTING. IT ALSO KNOWS THE CURRENT HELICOPTER 
LOCATION AND HEADING THRU COMMON /XFLOAT/. SUBROUTINE 
INIRDR WAS CALLED AT INITIALIZATION TIME TO CREATE 
A LIST OF TARGET DATA WHICH IS MAINTAINED IN A CIRCULAR 
LIST AND IS STORED IN COMMON /TARGET/. ALSO MAINTAINED 
IS A SET OF POINTERS INTO THIS TARGET LIST. TGTSTR POINTS 
THE FIRST TARGET BEHIND THE RADARSWEEP AND NOT VISABLE. 

TGTEND POINTS TO THE LAST TARGET WHICH IS AHEAD OF THE RADAR 
SWEEP AND AHEAD OF THE LOOKAHEAD AREA. A LOOKAHEAD AREA 
IS DEFINED AHE.AD OF THE SWEEP TO INSURE THAT ALL TARGETS 
ABOUT TO BE VISABLE TO THE RADAR HAVE BEEN UPDATED WITH 
REGARD TO DISTANCE .AND HEADING. THESE TARGETS ARE UPDATED 
BY PRIMARY UPDATE WHENEVER THEY HAVE ENTERED THE LOOKAHEAD 
AREA. IF ALL TARGETS IN THE LOOKAHEAD AREA HAVE BEEN UPDATED 
THEN AN ALTERNATE TARGET NOT IN THE RADAR SWEEP IS UPDATED. 
ONLY ONE TARGET IS UPDATED AT A TIME AND ONCE A TARGET IS 
VISABLE TO THE RADAR ITS POSITION WILL NOT BE UPDATED UNTIL 
THE RADAR SWEEP IS PAST IT. A TARGET UPDATE IS PERFORMED BY A 
CALL TO SUBROUTINE UPDATE. THIS CORRECTS THE HEADING AND 
DISTANCE INFORMATION AND INSURES THAT THE CIRCULAR LIST 
IS MAINTAINED IN THE PROPER ORDER. 

THE ALGORITHM STARTS BY CALCULATING THE ANTENNA HEADING 
AND THE HEADING DIFFERNCE OF THE TARGET POINTED TO BY TGTSTR. 
THIS INITIAL HEADING DIFFERENCE IS ALWAYS NEGATIVE AND 
HEADING DIFFERENCE IS CORRECTED FOR THE ANTENNA DIRECTION SO 
THAT A NEGATIVE HEADING DIFFERENCE MEANS THE TARGET IS BEHIND 
THE ANTENNA SWEEP. THE ALGORITHM PROCESSES EACH TARGET BY 
TRAVERSING THE CIRCULAR LIST IN ORDER OF ANTENNA DIRECTION. 
THIS PROCEEDS UNTIL A TARGET IS ENCOUNTERED PAST THE LOOKAHEAD 
AREA. EACH TARGET PROCESSED THAT IS VISABLE TO THE RADAR 
(IE. -TQTWDT < HDGDIF < TGTWDT ) GENERATES A NUMBER OF 
INTENSITY TRANSITIONS AT OFFSETS FROM THE TARGET DISTANCE. 
THESE CHANGES ARE STORED IN A LINKED LIST BY SUBROUTINE INSERT 
AFTER ALL TARGETS HAVE BEEN PROCESSED THE LIST IS TRAVERSED 
TO GENERATE A LIST OF INTENSITY CHANGES IN ORDER OF DISTANCE 
TO BE SENT TO THE RADAR HARDWARE. THIS OUTPUT DATA IS SENT BY 
A CALL TO SUBROUTINE SRDI. 

THIS ENDS THE PROCESSING DONE BY SUBROUTINE CREATE. 
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)EF INIT IONS: 

' INDIVIDUAL TARGET DATA 
NUNTGT - NUMBER OF TARGETS DEFINED 

TGTL0C<2, TGTID) - X AND Y LOCATION OF TARGET IN FEET 
FROM PRIMARY TARGET 

TGTHDG<TGTID> - TARGET HEADING RELATIVE TO TRUE NORTH 
IN RADIANS 

■TGTDSTiTGTID) - TARGET DISTANCE FROM HELICOPTER IN 
NAUTICAL MILES 

TGTPTR(2. TGTID) - INTEGER POINTERS TO THE NEXT 
COUNTERCLOCKWISE AND CLOCKWISE TARGETS 
TGTTYF(TGTID) - TYPE OF TARGET (0 = NON-BEACON, 1 = BEACON) 
TGTLGGtTGTID) - LOG 10 OF TARGET DISTANCE FROM HELICOPTER 
TGTTLT(TGTID) - ANGLE OF TILT FROM HELICOPTER TO TARGET 
IN R;ADIaNS 


: /STATUS/ STATUS OF POINTERS INTO CIRCULAR LIST OF TARGETS 

C TGTSTR - POINTER TO THE FIRST TARGET BEHIND THE RADAR SWEEP 

: AND NOT VI SABLE 

C TGTEND - POINTER TO THE FIRST TARGET AHEAD OF THE SWEEP 

C AND BEYOND THE LOOKAHEAD AREA 

C LSTDIR - LAST DIRECTION OF RADAR SWEEP (0=CCW, 1=CW) 

C NXTUPD - NEXT PRIMARY TARGET TO BE UPDATED. ALWAYS IN 

C THE LOOKAHEAD AREA OR EQUAL TO TGTEND 

C NXTALT - NEXT ALTERNATE TARGET TO BE UPDATED. ALWAYS 

C OUTSIDE OF THE VISABLE AND LOOKAHEAD AREAS. USED 

C WHENEVER NO PRIMARY TARGET IS READY FOR UPDATING 


C /3INTRE/ LINKED LIST OF INTENSITY CH.ANGES 

C BINLST( 1, NODE) - DISTANCE ALONG VECTOR OF THIS INTENSITY 

C CHANGE IN RADAR TARGET POSITIONS 

C BINL3T(2, NODE) - INTENSITY CHANGES OF LOWEST COLOR 

: BINLSTCS, NODE) - INTENSITY CHANGES OF MIDDLE COLOR 

C BINLST(4, NODE) - INTENSITY CHANGES OF HIGHEST COLOR 

C BINLST(.5, NODE) - POINTER TO NEXT NODE (LONGER DISTANCE) 

C HDBL3T - NODE CLOSEST IN DISTANCE 

C NXTNOD - NEXT NODE AVAILABLE FOR INSERT 

C /ANTENA/ ANTENNA PARAMETERS 

: ANTDIM - DIAMETER OF RADAR ANTENNA IN INCHES 

C RADOME - BEAM SPREADING FACTOR OF ANTENNA 

C ANGAIN - MAXIMUM GAIN OF ANTENNA IN DB 'S 

C 3RKPNT - DISTANCE AT WHICH STC CURVE CHANGES TO 

C 40 ^ LOG (DISTANCE) IN MILES 

■; ANTLOG - LOG 10 OF ANTDIM 

C /XFLOAT/ SIMULATION DATA 

C A(S) - PSIR - HELICOPTER HEADING IN RADIANS 

C A(103) - XPR - HELICOPTER X LOCATION RELATIVE TO 

C PRIMARY TARGET (IN FEET) 

C A(104) - YPR - HELICOPTER Y LOCATION RELATIVE TO 

C PRIMARY TARGET (IN FEET) 

C A(i05) - HPR - HELICOPTER ALTITUDE RELATIVE TO 

; PRIMARY TARGET (IN FEET) 

C A(35S> - D2R - DEGREE TO RADIAN CONVERSION 

C A (359) - R2D - RADIAN TO DEGREE CONVERSION 

C A(460) - FI - 3. 14159 

C A (461) - TWOPI — 2 PI 

C A (464) - VWIND - VELOCITY OF WIND IN KNOTS 

C A(46S) - PSWR - DIRECTION OF WIND IN RADIANS 
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PR 


u 


<?RAM 


tR I >^i: LE DEF I N I T I DNS 


AHDLEN 
ANTER P 


COMP-JiciD LtNGTH Or 
IN MILES 

PERCENT GF TARGET IN 
ERROR BETWEEN ANTENNA 


S5ET AHEAD Of- AC i UAL DISiANCE 

FRONT GF ACTUAL PG3ITIDN 
HEADING AND TARGET HEADING 


iiN L'C.vrH t!Z.a 

ANTGaN - ACTUAL ANTENNA GAIN IN DE CONSIDERING ANTENNA 
POINTING ERROR 

ANTHDG - HEADING OF ANTENNA IN RADIANS RELATIVE TO TRUE NORTH 
ANTING - CHANGE IN HEADING (IN RADIANS) FDR EACH 
ANTENNA POSITION INCREMENT 
ANTLOS - ANTENNA LOSS IN DB DUE TO ANTENNA ERROR 
ANTT3L(25> - TABLE OF ANTENNA LOSSES IN DB DUE TO ANGULAR ERRG 
ONE ENTRY FOR EACH DEGREE OF ERROR 
ASTLEN - COMPUTED LENGTH OF TARGET BEHIND ACTUAL DISTANCE IN 
MILES 

A3TPCT - PERCENT GF TARGET BEHIND ACTUAL TARGET POSITION 
ATK - ANTENNA TILT CONVERSION CONSTANT FROM RADAR 
UNITS TO DEGREES 

3CNGAN - BEACON GAIN SETTING IN DB. RANGES FROM -37 TO +3 DB. 


BECDEL — 
BGK 

CQEFl C2) 


C0EF2<3) 


BEACON TRANSMITTING DELAY IN MILES 
BEACON GAIN CONVERSION CONSTANT FROM RADAR 
UNITS TO VOLTS 

- CORRECTION FACTOR FOR TARGET LENGTH CONSIDERING 
PULSE WIDTH. 

(1) =2.35 USEC PULSE WIDTH 

(2) = .6 USEC PULSE WIDTH 

- PERCENT OF TARGET LENGTH FOR THE THREE INTENSITY 


LEVELS. 

(1) = PERCENT LENGTH OF LOWEST INTENSITY 

(2) = PERCENT LENGTH OF INTERMEDIATE INTENSITY 
•3) = PERCENT LENGTH OF HIG.HE3T INTENSITY 

C0NVRT(6) - ARRAY TO CONVERT DISTANCE (IN MILES) TO ENCODED 
RADAR DISTANCE FOR EACH RANGE SETTING 
CURTGT - TARGET ID OF CURRENT TARGET BEING PROCESSED 
DB - COMPUTED DE LEVEL FOR EACH INTENSITY LEVEL 
DELVL - COMPUTED RETURN SIGN.AL STRENGTH IN DBM 
DBZERO - TRIGGER LEVEL FOR SEA CLUTTER MODEL 
FIXFCT - CONVERTS FIXED SEA CLUTTER LEVEL TO VALUE USED TO 


COMPUTE MINRANGE AND MAXRANGE (I.E. CONTROLS 
CONTROLS CLUTTER LENGTH) 

FT2NM - FEET TO NAUTICAL MILE CONVERSION FACTOR 
HDGDIF - DIFFERENCE BETWEEN TARGET AND ANTENNA HEADING 
FOR CURRENT TARGET 

HDGFCT - FACTOR APPLIED TO SEA REFLECTIVITY FOR 
HEADING REALATIVE TO THE SEA 
I - WORK INDEX 
I ANT - ENCODED ANTENNA TILT 


lAWPCI - ENCODED ANALOG WORD RECEIVED FROM RADAR 
lAZM - ENCODED ANTENNA DIRECTION 
IBBR - ENCODED RADAR MODE 

1 = RADAR ONLY 

2 = BEACON ONLY 


3 = RADAR AND BEACON 


IBGN - ENCODED BEACON GAIN 

ICOMP - COMPLETE BIT VALUE 

IDAT.A( i-100) - ARRAY FOR OUTPUT DATA FOR 

IDIRF - CURRENT DIRECTION OF RADAR SWEEP 


RADAR HARDWARE 
(0=CCW, 1=CW) 


IDWPCI - ENCODED DISCRETE WORD RECEIVED FROM RADAR 


IMAPS - WEATHER OR MAP MODE 
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1 


— INTENSITY ADDED rOR 5TR0^4QEST TARGET 
4) - INTENSITY CHANGES FOR EACH LEVEL 

ALSO USED AS INTENSITY TOTALS FDR VECTOR CREATION 
4) - F^4CODED VALUE FOR ALL 4 INTENSITY LEVELS 

- CURRENT OUTPUT INTENSITY LEVEL 

- TOTAL INTENSITY LEVEL 
<0-l) = BLACK 
(2-3) = LOU'EST INTENSITY 
(3-6) = INTERMEDIATE INTENSITY 
(7- ) = HIGHEST INTENSITY 

PRESET RADAR INTENSITY IF = 1 

- PULSE WIDTH 

1 = 2. 35 US 

2 = .6 US 

ENCODED RADAR GAIN SETTING 

- CURRENT RANGE SETTING 

(1-6) = (2.5. 5.0. 10.0, 25.0. 

- RANGE AT WHICH INTENSITY LEVEL 

- RANGE AT WHICH INTENSITY LEVEL 

- SEED FOR RANDOM NUMBER GENERATOR 

- SEED FOR RANDOM NUMBER GENERATOR "7 

- SIZE OF LOOKAHEAD AREA (IN RADIANS FROM CURRENT 
ANTENNA HEADING) 

- INDEX INTO ANTENNA TABLE - 

- LAST INTENSITY LEVEL OUTPUT 

- MASK FOR ANTENNA TILT (ANALOG WORD) 

MAX - MAXIMUM ENCODED RANGE FOR CURRENT RANGE SETTING 
MAXINT - MAXIMUM INTENSITY LEVEL FOR CURRENT RADAR MODE 1 

MAXRNG(6) - MAXIMUM ENCODED RANGE FOR EACH RANGE SETTING I 

MEBR - MASK FOR RADAR/BOTH/BEACON MODE (DIGITAL WORD) 

MEGN - MASK FOR BEACON GAIN (ANALOG WORD) — 

MDIRF - MASK FOR ANTENNA DIRECTION (DIGITAL WORD) 

MDS(3, 2) - MINIMUM DETECTABLE SIGNAL IN DB FOR 3 INTENSITY ' 
LEVELS AND 2 PULSE WIDTHS _ 

MIN - MINIMUM ENCODED RANGE FOR CURRENT RANGE SETTING ; 

MINRNG(6) - MINIMUM ENCODED RANGE FOR EACH RANGE SETTING i 

MMAPS - MASK FOR MAP/WE.ATHER MODE (DIGITAL WORD) 

MFGR - MASK FOR PRESET GAIN (DIGITAL WORD) ” 

MRGN - MASK FOR RADAR GAIN (ANALOG WORD) i 

MRGSET - MASK FOR RANGE SETTING (DIGITAL WORD) 

M3ECT - MASK. FOR SECTOR TYPE 60/120 DEGS (DIGITAL WORD) 

NODE - CURRENT NODE FOR LINKED LIST TRAVERSAL 
NPT - NUMBER OF OUTPUT WORDS FOR SRDO 
NXTPGS - NEXT ANTENNA POSITION AZIMUTH 
NX.TTGT - NEXT TARGET IN DIRECTION OF RADAR SWEEP 
FROM CURRENT TARGET 

GHDDIF - HEADING DIFFERENCE OF LAST TARGET PROCESSED 

PCT - PERCENT OF TARGET LENGTH CONSIDERING DB LEVEL AND ' 

I NTENS I TY LEVEL ; 

PCTD3 - FACTOR MULTIPLIED BY DB TO INCRE.ASE TARGET LENGTH 
PCTMIN - MINIMUM TARGET LENGTH (AT 0 DB ) IN PERCENT 
PRSEA - RADAR ROWER RECEIVED FROM THE SEA 

PULSE - RADAR PULSE WIDTH IN USEC . j 

RANFIX - FACTOR FOR RANDOM EFFECT OF SEA CLUTTER 
RANGE! - RANGE AT WHICH INTENSITY LEVEL INCREASES (IN MILES) — 
RANGE2 - RANGE AT WHICH INTENSITY LEVEL DECREASES (IN MILES) ! 
RANGE5(6) - MAXIMUM RANGE (IN MILES) FOR EACH RANGE SETTING 
RANGR - CURRENT RANGE OF TARGET (OR SEA CLUTTER) IN MILES _ 
RDRGAN - RADAR GAIN IN DB (FROM -37 TO +3 ) 
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50. C. 100.0) MILES 
RISES (ENCODED) 
DROPS (ENCODED) 


I NT ADD 
INTCHG C 

INTENT ( 

INTOUT 

INTT07 


I PGR - 
IPULSE 


IRGN - 
IRGSET 

IRNGl 

IRNG2 

IRX 

IRY 

LCK.AHD 

LOOKUP 

LSTINT 

MANT 



RGK 


RHDO 

RNCCEF 

RNGLD5 

wC. AOii? 

SEAFCT 
SF AC T 


RADAR GAIN CONVERSION CONSTANT FROM RADAR 
UNITS TO VOLTS 

ANTENNA DIRECTION RELATIVE TO THE SEA 
(0 = HEAD SEA, 1 = DCWNSEA) 

TOTAL RANDOM EFFECT OF SEA CLUTTER 
LOG 10 OF CURRENT TARGET RANGE 
SEA RETURN POWER RELATIVE TO LEVEL 1 
THRESHOLD VALUE 

FACTOR WHICH MULTIPLIES WIND SPEED TO 
INCREASE SEA REFLECTIVITY IN HIGHER SEA STATES 
ANTENNA HEADING CONVERSION CONSTANT TO CHANGE 
RADAR UNITS TO DEGREES 


3IGSEA - RADAR REFLECTIVITY OF SEA PER UNIT AREA 
STC - STC LEVEL FOR CURRENT INTENSITY LEVEL 
3TCLVL - STC LEVEL OF HIGHEST INTENSITY SETTING 
STCORGiS) - OFFSETS ADDED TO STCLVL TO OBTAIN STC LEVEL FOR 
EACH INTENSITY LEVEL 

STHDDF - HEADING DIFFERENCE OF STARTING TARGET ID (TGTSTR) 
TGTLEN(6) - LENGTH OF TARGET FOR EACH RABGE SETTING (MILES) 
TGTWDT - HALF OF THE WIDTH OF TARGET (RADIANS) 

TILT - RADAR ANTENNA TILT (IN RADIANS) 

TLTDIF - TILT ERROR BETWEEN TARGET AND ANTENNA (IN RADIANS) 
TLTFCT - FACTOR MULTIPLIED BY TILT FOR SEA CLUTTER MODEL 
UPDTGT - TARGET TO BE UPDATED NEXT 

WNCGEF - FACTOR MUTIPLIED BY WIND FDR SEA CLUTTER MODEL 
WNDB - TOTAL FACTOR OF WIND, TILT AND WIND DIRECTION USED 
FOR SEA CLUTTER MODEL 
XANT - ANTENNA TILT IN DEGREES 
XBGN - BEACON GAIN IN VOLTS 
XRGN - RADAR GAIN IN VOLTS 
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''~323 • 

■ 324 i» 
S 2 ^ «i 


285 

2p6 

23 7 
233 
239 


291 

292 

293 

294 

295 

296 

297 
293 

299 

300 

301 

302 
30. 1 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 
313 

319 

320 

321 

322 

323 

324 

325 


C( 

000 

000 

000 

coo 

000 

000 

000 

con 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

QOC 

000 

000 

000 

000 

000 

oco 


AMERP = (ABS(W030IF) + A53 ( TLTD i f ) ) /D2'?^ 
AMfFR = ( aNTeRF^/RADSMe 4 ANrOIM/12,0) 
LB9<Up = Ipn<{ANTERR) + 1 

IF (Le3<UP *GT. 25) G9 T9 40 

A\:T|_9S 5 AN'TTaE ( 5 

ANTGAN = ANGAIN - Ain1TL9S 


C* * 4* > * 

c*** 
c* 

c * » * 
c* 

c* * 

C* 

r ♦ 

C** * 

p4# ♦ 

C*** 

C** * 

C 
C 



RETURN P9WER (D0) = TRANSMIT PQwER - 
3Q4L9G10(4*PI ) - 40*L9310(RANGE IN METcrS) + 24(AMTE''JNA 
GAIN) + 204L9GlC(LA^iBDA) + IO^lSGIO ( S IGMA ) 

- T/R L3SS..f 

whfre; 

transmit PR.-jEr = 7000 WATTS » 6^.45 dEc^SElS*** 

antenna gain - 6.4 + 20Lv9G10HHHHH»L9Gl0( ANTOIM) - n9INTIN 
LAMriDA - o*03l9 METERS* •« L9S3 

sigma = 1000 squarf Meters Reelectivity. , . 

T/R L^SS = TvyO-wAy (USER INdjT)... 

THIS reduces T9: 

OBUVL = -95.13 - 40*L931Q(RANGE IN NAUTICAL MILES) 

+ 2* ( antenna GAIN) , , , 


IF ( ISBR .EG. 2) 39 T0 25 

03LVL = -95.13 + 2 .O 4 ANT 3 AN - 40*0*RNgLS3 - TRL9SS 
STCLVL = -85,40 + ANGAIN - 23.04RNGL9G - RDRGAN 

Ip (RaNGR .GE. RN3BRK) 

♦ STCLVL - -87.30 + 2 . 0000 * aNGA T N - 39.364RnGL9G - RDRGAM 


Ip ( I M ApS 


■ • 


1 ) 
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32ft 

327 

328 

iig 

330 

331 
333 

333 

334 

335 

336 

337 

338 

339 

340 

341 
34? 
34 3 

344 

345 

346 

347 

348 
34<? 
35u 
35l 
35p 
353 
35^ 

355 

356 

357 
258 

359 

360 
36 1 
362 
36 3 

364 

365 

366 


326. 

327. 
323. 
329, 
33Q. 
331 

332. 

333. 

334. 

335. 

336. 

337. 
333, 

339. 

340. 
341 

342. 

343. 

344. 
343 , 

346. 

347. 

345. 
34g, 
330- 
351 
352' 

353 . 

354. 

355. 

356. 

357. 
353 . 

359 ' 

360. 

361. 

362. 

363. 

364. 

365 . 

366. 


000 

OUO 

000 

000 

000 

000 

000 

ooc 

000 

000 

ooc 

000 

000 

000 

ooc 

000 

000 

000 

000 

COC 

000 

coo 

000 

000 

000 

000 

000 

300 

coo 

000 

000 

000 

000 

000 

000 

ooc 

000 

ono 

000 

000 

000 


♦ STCLVL e -92,33 4 1 ,3333*ANiaAIN - 29,gO*R\aLeQ - RDRGAnj 

MAXIMT = 3 
IF ( IB69 ,rQ, 3) MAxHt * 2 


08 25> j » li'^AylMT 

If" (SjC ,Lr. '"DS( J, iFJLSfr) ) 
IF ( DB ,Lr. 0.0) 33 T6 25 


STc = stclvl + srcaRGfj) 
StC « M0s( J, IPULSE) 

OB = . 0 =lVl ” 3 rc 


PCT = (PCTDB*D3 + PCTM IN ) *C9EF2 U ) 
Ra^^GrI » RANQR - PCt^aHOLrN 

IRNqI = IFIX{RANrF1*C»NVRT( IR^GET) ) 

IF { IRNGl .GE. MAX) 39 T9 25 

IF (IRNGl ,LT. MIN) jRNGI = Mjkj 

RAN3E2 = ?CT»A3TLEN + HANQR 
I.PNG2 s IFIX(RANGE2*C0 NvRT( IRGSeT) ) 

IP ( IR^'G2 ‘LT. MIN) q9 T9 25 

IF ( IRNG2 .GE. .MAX) IRNG? = MAX - 1 

InTcPQ(2) = INTcH 3(3) = lNrcHG(4) = 0 
INTCHG(J) =-i 
InTCHG(J4i) » 1 

call INSFRT( 1RNg 1> IRN,32i INTCHy) 

25 C9NTINUE 


c 

c** ♦ 
c** * 


n 


n 




BEa^SN TArGETG... 

RETURN P0.vER (Db) = TRANSMIT P'3WER - 
HO^LeGlo ( 4*PI ) - 20*LeGlo IPANG e IN meTERS) + 

ANteN'NA GAI'N 4 204LBQ10(LaMB0A ) - SN^.wAY L9S3 
4 bEacsn antenna gain... 

where ; 


n 


n 


98 



367 


367 

•000 

36A 

• 

363 

•000 

36^ 


36g 

,000 

370 

m 

370 

• 000 

371 

«* 

371 

• OGC 

3?2 

m 

372 

.000 

373 

m 

373 

• 000 

374 

- 

374 

• oon 

375 

- 

375 

•OQO 

37o 


376 

• 000 

377 

• 

377 

• GOG 

378 

m 

373 

,000 

373 

m 

379 

•000 

38q 

m 

330 

•000 

331 


381 

.000 

382 


332 

• 000 

383 

W 

33 3 

• 000 

384 

m 

334 

•000 

383 

- 

385 

• 000 

386 


336 

• 000 

337 

- 

387 

• 000 

388 

m 

333 

• 000 

383 

m 

389 

• 000 

330 

- 

390 

.000 

391 


391 

• oon 

392 

m 

392 

•000 

3g3 

m 

3g3 

.000 

394 

m 

394 

• 000 

395 

« 

395 

.000 

396 

m 

396 

• 00 c 

397 

m 

397 

.000 

298 

m 

298 

,000 

399 

m 

399 

• Ooc 

400 

- 

400 

.000 

401 


401 

•000 

402 


402 

.000 

403 

• 

403 

.000 

404 

w 

404 

• 000 

405 

w 

405 

•000 

406 

n 

406 

•000 

407 

m 

407 

.000 


c ^ 

transmit power 

c* 

AN'tfNiNA gain 

c* * * 

L A ^ B 0 A 


■one-way L0SS 

c** * 

RFACO'N aNtfNNA gain 


THIS REDUCES TO: 


OrLVL = -56.16 

c * 

StCLvL = ,51,00 

c * 

00 = DrLVL 


DB ~ “5»l6 




) , PSlNjINr, L33S 
RETUR.Ni pULrE 
p-REQUENCy) • • • 


C 

C 

c 


c 

c 

c 

c 


+ ANTE.nINA gain - 0.5*(T/R L8SS))..t 


■ ***■*»*#* 


IF ((IflBP .£3. 1) *cR. (TGTTYP(CURTGT) .EQ. 0)) 

Gc T3 40 

C8 = -5,16 + ANTGAN - 0*5»TRL9S3 + BCNGAN 

IF (D3 ,LT. OtO) G6 T3 40 

RANqEI = RANqR + bECDEL ahclen 
IRNGI = IFIX(RANGE1*C9NVRT(IRG3ET) ) 

If (IR^'Gl .GE. mAX) G3 T8 40 
IF {IRNGI ,Lr. ^IN) IRNGl - MIN 

RANGES a RANGR + BECDEL + aSTLEN 
IRNG2 = lFIx(PANGE2*CaMVRT( IPGS eT) ) 

IF (IRNG2 ,Lr, MIm) G0 T0 UQ 
IF (IRNG2 ,GE. MAX) IRNGE = MAX - l 

InTcH6 (2 ) = INTcHgO) = 0 
I\'TCHG(4) = 1 


CALL INSERT( iRNGli IRNG2^ InITCHG) 



’ 1 


408 

m 

408 

• ooc 

C ** * 

target is In tme lsok-aheao area... 


409 

w 

409 

• 000 

c 




41.0 

m 

410 

• Ooc 


35 

IF (T3TEND ‘Eo. CURTGT) TGTEND = NxTTgT 


^11 

m 

411 

•ooc 



IF((NjxTL'PD »NE. CURTGT) .or. CJRdTGT .nE* 0)) 


412 

m 

412 

,ooc 


* 

G6 jO 4O 


413 

m 

413 

.000 

c 




41 4 

m 

414 

• 000 



UPDTGT » CURTGT 




415 

• oco 



nxtupd = nxttgt 


416 

« 

416 

• 000 

c 




417 

m 

417 

• 000 


GET THE NEXT TARGET... 

i 

418 


413 

.000 

c 




419 


419 

• 000 


40 

CUR'^G^ ® NiXTTnT 

n 

1 

42n 

mm 

420 

•000 



G6 T9 10 

42i 

m 

421 

•000 

c 



- j 

422 

m 

422 

•000 


target is In FReNT OF THE L69<-AHEAD AREa... 


423 


423 

.000 

c 



i 

424 

m 

424 

• 000 


45 

lSTdIR * tOlPF 

. j 

425 

m 

425 

• 000 

c 




426 

m 

426 

• 000 



IF (CURTGT *Eg. TGTEND) G9 T9 90 

1 

1 

427 

m 

427 

.000 



IF (NXTUPD ,EC. TGTEND) NXTUPO 3 CURTGT 

428 

m 

423 

•Ooc 

c 



} 

429 

m 

429 

• Ooc 



TGTEnO = CURTGT 


430 

W 

430 

• 000 

c 




431 


431 

• ooc 

c^ 


find an AI.TF.RNATE TARGET... 

. 1 

432 


432 

• 000 

c 




433 

• 

433 

.000 


50 

IF (UPDTGT .NE. 0) G9 T6 55 

"1 

1 

434 

mm 

434 

,000 

c 



435 

m 

435 

.000 



Ic((TgTENC .EC. TGTPTRdDIRc + 1. TGTSTR)) .AND. 

1 

436 

• 

436 

.000 


* 

{STHDDF .LE. L9KAHD - TWRPI)) G9 T9 57 

n 

T 

437 

m 

437 

• Ooc 

c 


438 

W 

438 

• 000 



yPOTGT = nXTAi_T 

1 

439 

m 

439 

*000 



NXTAlT 3 TGTPTR(2i UPDTGTj 


44o 

W 

440 

•000 

c 




441 


441 

• ooc 



if((idirf ,eo. 1) .And. (updtgt .eq. tgtstr)) 

1 

1 

442 

m 

442 

,00c 


* 

NXTALT = tGtfnD 


443 

w 

443 

• 000 



I|7((lDlRF .Eg. 0) .And. (upDTgT .fq. T^TEmD)) 

^ 

444 


444 

•000 


* 

NXTALT = tgtstr 

) 

1 

445 

m 

445 

.000 



IF (UPOtGT .eq. TGTSTR) TGtSTR = TGTPTR(2 - IDIRF, 

UPnTG. ) 

446 


446 

• 000 



IF (UPOTGT .eq. TGTEND) TGTEND » TGT?TR(1 + IDiRFj 

UPOTGT ) 

447 


447 

.000 

c 




448 

w 

448 

• 000 


ijPdATE the TARGET... 



100 



i 

«* 

44g 

.OOC 

c 




■t 

450 

• 000 


55 

TGTIq » 

**51 

- 

451 

• OOC 

r 

w 


,■ 

r ^5? 

« 

452 

• 000 


CALI update ^ 

' 453 

* _ 

- 

453 

*000 

c 



454 

- 

454 

• OOC 


r* 453 

- 

455 

.000 



i 45^, 

w 

456 

♦ 000 


StA CLU''"rE:R m9qE|_ (IF 'TIlT' is GrEATFo T^AN + 2.0 

457 

m 

457 

• 000 

r 

w 


Degrees, ns sea clutter is oisplayeO) • . , 

r-453 

w 

438 

•000 

u 



, 459 

«■ 

459 

♦000 


57 

IR ((TILT .ar. 0•0^^) •OR^ (IBBR .EQ. 2n GS T0 6 q 

463 


460 

• 000 

c 



461 

m 

461 

♦ 000 



TILTO = TILT*R2D 

^ 463 

m 

462 

• 000 



IF (TILTO .GT. -1.0) TILTD s -1.0 

; 46 3 

m 

463 

• OOC 

c 



464 

m 

464 

•000 



RANGR ' (-HPR*FT2^M)/(TILTD*DHR) 

,—465 

w 

465 

,ooc 



RNGLftG = Lc)GIQ(RaNGR) 

■ 4 66 

m 

466 

•000 



STCLvL = - 85,4 + AN 3 AIN - 23 * 0 »RMQLe 3 - RDR 3 AN 

467 

m 

457 

•000 

c 



46^? 

m 

455 

•000 



IR (RaNGR .GE, RN3BRO 

: 469 

m 

469 

• 000 



3tCLVL = .87.80 + 2.0000-» aNGa I^ - 39 . S6*RNGL9G 

■ 470 

m 

470 

• 000 


* 

- ROR-jA'M 

471 

m 

471 

• 000 



IR (IMAPS ,£3, 1) 

—472 

m 

472 

,000 



StCLvL = .y2,33 + 1 . 3333 *aNGA I N , 2g,gO*RNGL9G 

^ 473 

- 

473 

• 000 



« RDR'GAN 

' 474 

• 

474 

• 000 



STC = STCLVL + STC5RG ( 1 ) 

475 

w 

475 

,000 



IR (SjC ,Lt. M0S(1, IPJLSe)) STC = MDS( 1^ IPULSe) 

r 476 

w 

476 

• 000 

c 



' 477 


477 

• 000 



pulse = 2.35 

475 


478 

• 000 



IF (IPULSE .Eo. 2) PUlSE = 0 . 6 Q 

'~4 79 


479 

•000 

c 



: 430 

- 

480 

• ooo 



IR (ANTHDQ ,GT. PI) ANT^DG = ANTHDG - T.mOPI 

431 

- 

481 

.000 



PSW^E^ = PS'*^R 

,—4-32 

m 

482 

• OOC 



IR (pSi^NEw •QT.pI) pS^iNiEw = pSwMEw " T,,^,9pl 

: 433 

m 

483 

•OOC 

c 



' 484 

- 

434 

•000 



RHOG = ABS( ANTHDG - PSWNEW)/PI 

435 

- 

^85 

,000 



IR (VUInO ,Lr. 5,0) RhDG = 0,0 

;~436 


486 

• 000 

c 



: .437 

«■ 

437 

• 000 



SIqSEA = -40.00 - TLTFCT*TILTD - HDgFCT*RHDg + 3EAFcT»v'^IN 

433 

» 

488 

• 000 



pwRSEA = -66.72 + 2.Q*ANGAIN + 10 »0* ( LSGlO ( PULSE 1 ) 

—438 

m 

439 

,000 


•» 

30,0*RNGL9G - TRLOsS - 10,0*ANTLoG 
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- 

490 

•000 


+ SI3SEA 


491 

* 

491 

.000 


3EADB = PAt^SEA - 3TC 


49P 


492 

• OCG 

c 


1 

49:^ 

m 

493 

• 000 


IP (SEADB .LE. O.Q) GB T5 60 

1 

494 

- 

494 

,000 

c 


49y 

w 

495 

• 00c 


= SEAD3*pIXpcT 


49<^ 

- 

496 

• 000 


pnCbef r (Panfix -*• ra'JDS*wnD3 ) *rangr 


497 

- 

497 

.000 

c 



*♦9^ 

- 

493 

.QQO 

c*** 

GET A SAmPl 5 3F uNIFe^MLY oISTrIBjTEc W^l T f NOISE... 


499 

m 

499 

.COO 

c 

1 

1 

500 

tm 

500 

.000 


Call randum( irx/ iRYy zG "J\ ) 

« I 

501 


50 1 

.000 

c 


502 

- 

502 

.000 


RANgEI = RA\'gR/<1.C + wNDB) + <ZGW"- - r,.5)*RNcOEF 


5C3 

m 

503 

.000 

c 


. ) 

5C4 

- 

504 

• 000 


GET A SAyPLE of UNIFdRMLY DISTRIBUTED WmtTe NBISE... 


5Q5 

m 

5 O5 

.000 

c 


n 

506 


506 

.000 


call rANDumI IrX< IrY/ 7-G.Y'-) 

1 

507 

m. 

507 

.000 

c 



503 

m 

503 

•000 


RANGE2 = RAnGR*(1.0 + ,\NDB) + (ZGWN - O.5)"*RNC0EF 



50 q 

- 

509 

.000 

c 


1 

510 

m 

510 

.000 


If (RANgEI .gT. RANGES) QB TO &0 

1 

511 

- 

511 

.000 

c 



512 

- 

512 

• 000 


IRNgI = RAngF1*cB^JVRT( IRgSET ) 

n 

513 

- 

513 

.000 


IRNG2 = RA'gGE2*C0NvRT( IRGSET ) 

1 

5l«. 

- 

5 I 4 

.000 

c 



515 

w 

515 

.000 


If (IRNGI .GE. mAX) G0 to 60 


516 

m 

516 

.000 


IF (IRNGI .Lr. '^IN) IRNGI - YIN 

1 

517 

m 

51 7 

•000 


IF ( IRNG2 ,LT. MIN) Ge) T9 60 

■ 

513 

m 

513 

,000 


IF (IRNG2 .Ge. max) IRnG2 = MAX - 1 


519 

m 

519 

•000 

c 


n 

520 

- 

520 

• 000 


INT0HG(2) = 1 

j 

521 

** 

521 

.000 


INTCHGO) = 0 


522 

m 

522 

.000 


INTcHG(4) = 0 


523 


523 

.000 

c 


i 

524 

■t 

524 

• 000 


CALI. INSERT( IRNGI > IRNG2. InTCHG) 


525 

m 

525 

.000 

c 


n 

526 

- 

526 

• 000 

c * * »*•»■ 

create the vECTBR... 

527 


527 

• coo 

c 


j 

523 


523 

• 000 

60 

N9DE = HDBL3T 


= 29 

m 

529 

.000 


LSTINT = NPT = 1 

n 

530 

m 

530 

.000 


IDATA(l) = INTENT(LSTInT) 

\ 

i 


~j 
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531 

• 

531.000 

532 

- 

532.000 

533 


233*000 

534 

- 

534.000 



5 3 5 • 0 0 

536 

- 

536 . OOC 

537 


537. 000 

536 

• 

-33.QG0 

53 g 

m 

53.3,000 

54 n 

m 

540.000 

541 

m 

541.000 

542 

m 

542,000 

543 

m 

543.00c 

244 

m 

544.000 

545 

m 

545 * 000 

546 

•t 

546,000 

547 

m 

547.000 

546 

- 

548.000 

54-4 


543.00c 

350 

- 

550 .000 

551 

- 

551 .000 

552 

m 

552 . 000 

553 

m 

553,000 

554 

- 

554.000 

555 

- 

555 . 000 

556 

• 

556.000 

557 

m 

557,000 

556 

m 

558 * 000 

559 

m 

559*000 

560 


560 ,000 

561 

- 

561 » 00 C 

562 


562 *000 

563 


563.000 

064 

m 

564 . 000 

565 

- 

565.000 

066 

m 

566.000 

567 


567,000 

566 

m 

563.000 

56 S 


569.000 

57 o 

- 

570.000 

571 

- 

571.000 

572 

• 

572.000 

573 

m 

573 * 000 

574 

•* 

574.000 

575 


575 • 000 

576 

- 

576.000 


c 

C' 

c 


c 

c 

c 


IMCHG(I) = IMCHG(2) = 
Gfl Ta 75 

create thp: vector... 


I\'TCHG(3) = 0 


70 ' I\TchG(1) = IN'TchG(I) BlNi|_ST(2/NeQE) 

I\'TCHG(2) . INTCH3(2) ^ B I\LST ( 3» ^3DE ) • 

I\TCHG(3) = r\TCHG(3) + B I NL ST ( 4 / NODE ) 

iNijTaT = INtCHG(I) + 2*I\'tCHG(2) + 3 * I \'tCHG ( 3 ) 

IN.TADD = .1 

IE (INTCHG(2) ,\E, 0) INTAOD = 2 
Ir (I.VTCHg(3) .NE. 0). INTAOD = 4 

INTT3T = I'MTTST + IMTADD 

INtoUt = I^iTTST/2 + 1 

IE (IN'TTST «EQ. 6, INT3UT = 3 

I..AXINT = 4 

IE ((I3BB .EC. 3) .and. (INTCHG(3) .EC. 0)) MAXINT = 3 
IE (iNTSUy ,Gt, MAXINJ) INjOUT » MAXlMj 
Ip (INTOUT .EC. LSTINT) G0 T0 74 

NPT = NPT + 1 

IOAtAINot) = If3R( INteNTC INT9UT).8INLST( ) 

LSTINT = IMT0UT 

N0DE = 3INLST(5#N0DE) 

Itr (N0DE .NE. 0) 30 T9 70 

CREATE the C0MPLETE VECTOR AND OUTPUT THr DATA... 


- 74 
C 

75 

C 

c** * »■» 

c 


c 

c 


NPT = NPT f 1 

IDATA(NPT) = I9R(IC0MP. MAXRNG ( I RGSET ) ) 
CALL SRD0( IDaTAjMpT) 

999 RETURN 


C* * » * ♦ » 

c 

ENq 


END OE subroutine 'LREATE'... ****** 
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SUBROUTINE INSERT 
DESCRIPTION: 

INSERTS INTENSITY CHANGES INTO LINKED LIST 
LIST ARRANGED IN ORDER OF RANGE. THE SUBROUTINE 
IS PASSED AN ARRAY OF INTENSITY CHANGES FOR LEVELS 
0, 1, 2, and 3 RESPECTIVELY. ALSO PASSED ARE 2 
RANGES. THE FIRST IS WHERE THE INTENSITY CHANGES 
OCCUR WHILE THE SECOND IS WHERE THE INTENSITY 
LEVEL REVERTS BACK. FOR A CHANGE OF LEVEL 1 TO 
LEVEL 2, THE ARRAY WOULD CONTAIN A -1 FOR LEVEL ^ 
AND A +1 FOR LEVEL 2. THIS INDICATES 1 LESS TARGET 
AT LEVEL 1 AND I MORE AT LEVEL 2. IF INTENSITY 
LEVELS ARE INSERTED AT A PREVIOUSLY INSERTED RANGE, 
THE INTENSITY LEVELS ARE ADDED. OTHERWISE, A NEW 
MODE IS ADDED. 


COMMON DESCRIPTION: 

SEE SUBROUTINE CREATE 


VARIABLE DESCRIPTION: 

JHEAD - CURRENT NODE IN LINKED LIST 
INSNEG - SET TO 1 AFTER NEGATIVE INTENSITY CHANGES HAVE 
BEEN PROCESSED 

INTCH(4) - INTENSITY CHANGES BEING INSERTED 

INTENC(4) - LOCAL ARRAY TO HOLD INTENSITY CHANGES 

IRNG - CURRENT ENCODED PJ\NGE BEING INSERTED 

IRNGl - ENCODED RANGE FOR INTENSITY CHANGES 

IRNG2 - ENCODED RANGE FOR REMOVAL OF INTENSITY CHANGES 

NEWNOO . INSERTED NODE 

OLDHD - PREVIOUS NODE IN BINARY TREE 


SUBROUTINE CALLED: 


NONE. 



I ‘ 


namf: 

1 

3 

4 

n 

6 

7 

9 

10 

11 

12 

13 

14 

15 

U 

17 

18 
19 
ao 
ai 
22 

23 

24 

25 

26 
27 
2a 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


S; INSERT. HFLRAD 


1 

2 

3 

4 

5 

6 
7 
3 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 
,37 

38 


000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

ooc 

000 

000 

000 

000 

ooc 

000 

000 

000 

000 

000 

000 

ooc 

000 

000 


c»-» *♦* 

* 

C ** *•) 

* i 

c 

c*** 

c*»* 

c 
c< 
c 

c 
c^ 
c 

c 


T I TLE • . • 
LAST modified,: 


s:insert^ helrad 

20 JAN^ *32; VOBEN 
• » • 


»*■»*■»* 
♦ •*»*** 


this subroutine inserts intensity (COLOR) changes into 

THfT binary jReE arranged BY RANGE VALUES... 


SU8R0UTINES called: NONE*.. 


^♦*»*»*-* 

c 


Subroutine insertc Irngi/ Irnge. iNTCfiG) 

*♦*♦»*#**»♦■**•»**• 

COMMON BLOCKS* . . 

.■»»♦*♦** *#*!.*» 

common 


#****•» 


* 

* 

♦ 


/BlNTRE/aiNLST(5il50)/ HDBLSTi NXTNQd 
/TARGET/TGtL0C(2^ 100), TGTDST(IOO), tGtHDG(IOO), 
TaTpTR{2, 100)i TgTLOG(IOO)^ TGTTLT(IOO). 

TGtTYP(IOO), 

NbRTGT/ NUMTGT/ tgtid 

/STATUS/LSTOIR# NXTAlT^ NXTURDj tgtend, tgtstr^ LaZM 


c 


integer 


integer 


declarations. . • 

BiNLSTi HDBLSTi OLDHD/ TGTTYP 
lNTCHG(4)i INTENc(4) 


*»***■» 
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) 


3 9 - 

33 

000 

40 • 

40 

000 

41, - 

41 

000 

42 - 

42 

000 

43 - 

43 

000 

44 " 

44 

000 

45 - 

45 

000 

46 - 

46 

000 

47 - 

47 

000 

48 - 

48 

000 

49 - 

49 

000 

5q • 

5o 

000 

5l - 

51 

00c 

5? - 

52 

000 

S3 - 

53 

00c 

54 •* 

,5^ 

000 

55 ■ 

55 

000 

56 " 

56 

00c 

57 - 

57 

00c 

58 - 

58 

00c 

59 ” 

59 

000 

60 • 

60 

000 

6 1 " 

61 

000 

62 • 

62 

000 

63 " 

63 

000 

64 - 

64 

000 

65 • 

65 

000 

66 *• 

66 

000 

67 ' 

67 

000 

63 ^ 

68 

000 

69 - 

69 

000 

7o - 

70 

000 

71 - 

71 

000 

7? - 

72 

000 

73 - 

73 

000 

74 - 

74 

000 

73 • 

75 

000 

76 w 

76 

000 

77 ■. 

77 

000 

78 - 

73 

000 

7g - 

79 

000 


c ^ 

c . _ 

c*»* INITIALIZE and CeNVERT THE RANGE* *t 1 

C J 

SLDHO ’ D 

JHEAO a HDBLST 

IRNG = IRNGi ! 


INTENC(a) » INTCHG(2) 
INTE,\C( 3 ) » INTCHGO) 
InTEnC(A) s INTCHGIA) 
INSNEG a 0 


C*** CHECK NeDE»*. 

c 




10 If (Jheaq .eq. 0) Ga la 50 

If ( IRNG *- 8INI.ST( 1/ jHEaD) ) 50^ 30> 40 

c 

c*^* same range - AOO intensity... 

c 

30 BINLST(2^ JhEAq) * 9 1 NlST ( 2i J hEAD ) + INTENc(2) 

BINLST(3< JHEAO) » B I NLST ( 3/ JHEAD ) + INTENC(3) 
3INLST(4J JHEAD) = aiNLST( 4 i JHEAD) + INTENC( 4 ) 
GB TB 60 
C 

c»** get the next point... 

c 

^ldho = jhEao 

JHEAD a BINLST{ 5^ BLDHD) ; G8 TO 10 

c 

C*** insert a new N0DE... 

c 


1 

. 1 


1 


) 




~l 


"1 


50 


NEWNBO = NXTN80 
NXTNBO a NXTNBO 
BINLSTC 1/NEWNBD) 
RlNLSTOiNEWNSQ) 
giNLST(3iNEWN0O) 
SINLST(4#NEWN6D) 
aiNLSTI S jNEWNBD) 


+ 1 

a IRNG 
a INTENC(2) 
a INTENC(3) 
= INTENC(4) 
» jhead 


“1 


I 

i 

1 


1 

1 



i 
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SO 

f* 

80 

• 000 

c 


31 


31 

• 000 

c* 

♦ ♦ IJO^ATE the POiNTEpS... 

82 

m 

32 

• 000 

c 


S3 

«• 

33 

•000 


IP (BLDHD .EQ. 0) GO TO 55 

34 


34 

• 000 

c 

85 

m 

35 

•000 


BINLST(5/0LDH0) => NEWN80; GO TO 56 

36 

- 

86 

•000 

c 


37 


37 

•000 


55 HDBLST = NEvn'NOO 

3S 

m 

38 

.000 


56 JHEaD = NE'^iiNlflO 

89 


39 

•000 

c 


9o 

m 

90 

• 00c 

C*** INSPRT negative P9INT... 

91 

m 

91 

.000 

c 


92 


92 

• 000 


60 Ip (INSNEG .NE. 0) GO TO 999 

93 

m 

93 

• 000 

c 


9^ 

- 

94 

.000 


INTENC(2) * -INTENC{2) 

95 

m 

95 

.000 


INTE‘'C(3) = ,lNTrNC(3) 

96 

- 

96 

.000 


INTEMc(A) 3 -IMTENc(4) 

97 

- 

97 

.000 


INSNEG 3 1 

95 

m 

98 

•00c 


IRNG = IRNG2; GO TO 10 

99 

m 

99 

.000 

c 


ICO 

m 

100 

• 000 


999 rETurN 

101 

m 

101 

• 000 

c 


1C2 

" 

102 

.000 


103 

• 

103 

•000 

c#***# END 0F SUBf^eUTINE MNSERT*,.. 

104 

w 

104 

•OQC 


IO5 

m 

IO5 

.000 

c 


106 

m 

106 

• 000 


ENq 
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subroutine update 


DESCRIPTION; 

UPDATE UPDATES ONE TARGET'S DISTANCE AND SEARING 
AND rtOVES IT TO THE CORRECT LOCATION IN THE CIRCULAR 
TARGET list 

COMMON DESCRIPTION: 

SEE SUBROUTINE CREATE. 


variable DESCRIPTION: 

CUPTGT - CURRENT TARGET FOR CIRCULAR LIST TRAVERSAL 

HDGCHG - CHANGE r^4 HEADING FROM LAST UPDATE 

IDCCW - target id of COUNTERCLOCKWISE TARGET 

IDCW - TARGET ID OF CLOCKWISE TARGET 

LENTRV - NUMBER OF TARGETS TRAVERSED 

NEWHDG - NEW HEADING OF TARGET 

NXTHOG - HEADING OF CURRENT TARGET BEING TRAVERSED 
QLDHDG - HEADING OF TARGET LAST UPDATE 
TCTID - ID OF TARGET BEING UPDATED 


SUBROUTINES CALLED: 

HDGDST - COMPUTES HEADING AND DISTANCE OF TARGET 



NAMP 


StSCIUPO.wELRAD 


C* # -4 4*^ 

c*^-»** 


titlf:,,. 
last MtjDIFlED: 




c*** this Subroutine updates 8ne target's distance and 

c*** reaching A'Nq viavES IT T9 THE c^RR^CT LOCATION IN THE 
C*** CIRCULAR target LlST».» 
c* * ♦ 


S’SCIUPDj HELRAO **#h.** 

20 JAN^ *321 VOBEn «*#♦*** 


SUBROUTINES CALLED: 

HOqDST - C9MPUTES THE HEADINq AND DISTANCE OF EACH 
target. . . 


SUBR9UTINE update 


CSmmON 3L9CkS... 


♦♦**** 
• »*♦*** 


CSHweN 

► /XFLBAT/ A(500)/IFIXED/IA(25U) 

► /HRCeM / RH(50), IH(50)/KC(25) 
common 

^ /BINTRE/olNLST(5/l50)i HOBLST^ NXTNgo 

^ /TArgET/T3Tl5C(2^100)/ TgTdST(100)> TGThDG(100)> 

► TGTPTR( 100) j TGTLOG(IOO)/ TGTTLT(l0O)i 

' ' TGTTYP(100)i 

^ N3rTGT> NUMTgT/ TGTIQ 

» /STATUS/LSTDIR# NxTALT* vXTUPD, TGTE.'iD, TUTSTR, La^M 




39 


39 

• 000 


equivalences •• • ****** 


^0 

■I* 

40 

. 000 , 




m 

41 

.000 

c 



W?. 

m 

42 

.000 


EQUIVALENCE 


^3 

m 

43 

.000 

' » 

( A( 358 )j ( A( 460 )> PI)^ ( A'(<^ 6 l)< Tw0Pl) 

n 


m 

44 

.000 

c 



^5 

m 

45 

.000 


1 

_ « 

46 

m 

46 

.000 


DeCLaRATI9NS«*« ****** 


47 

m 

47 

.000 


n 

1 

43 


48 

.000 

u 


. i 

49 

m 

49 

.000 


SEAL 


So 

m 

50 

.000 


NEwHDGj NxTHDG 

n 

5l 

m 

5l 

.000 


integer 

1 

i 

5P 

m 

52 

.000 


CNTTGTi TGTENq. TGTIq/ FGTPTr, TGTSTqi TGTTyP 

- j 

53 

m 

53 

.000 

c 



54 

m 

54 

.000 


n 

55 


55 

.000 

c****^* ^ *^ *^* * *^^ ******** *^***^^^^*^*^****^*^*1h¥^ ^ it ******^1i^*1Hi.*^k 

) 

56 

m 

56 

.000 

c 



57 

m 

57 

.000 

c * ■* * 

CaMPUTE distance and bearing... 

■n 

58 


58 

.000 

c 


1 

1 

59 

m 

59 

.000 


SlDHOG = TGThdG(TGTId) 

J 

60 

•P 

60 

.000 

c 



61 

w 

61 

.000 


Call hogost 

1 

62 

« 

62 

• 000 


NEwHOG » TGTHDG(TGTID) 

. 1 

63 

w 

63 

.000 


HDGcHG f NEwHDG - 3i_DH0G 


64 

m 

64 

.000 


IE (HDGCHG .GT. PI) H03CHG s HDGCHG - TW 8 PI 

j 

65 

m 

65 

.000 


IF (HDGCHG .lT. -PI) HDGCHG a HDGCHG + j'MQPl 

. 1 

66 


66 

.000 


NE>(Hr)G a Qco^D^ ■*■ hdGCHG 


67 

«• 

67 

♦ OQO 

c 



63 


68 

.000 

c** * 

CHECK FeR PRaPES PLACE IN LIST... 

1 

69 


69 

.000 

c 


J 

70 

m 

70 

.000 


If (HD^CHG *lT. 0*0) go to 40 


71 

« 

71 

• 000 

c 


» 

72 

m 

72 

.000 

c#** 

traverse CL 8 C<WISE... 

1 

. 1 

73 

w 

73 

*000 

c 



74 

m 

74 

• 000 


CNTTGT a tgtio 

n 

75 

m 

75 

.000 


LENTRV - 0 

1 

1 

76 

m 

76 

.000 

c 


- 

77 

m 

77 

.000 

30 

LENTRV a LENTRV + 1 


78 


73 

.000 


CnTTGT a TGTPTR(2iCNTTGT) 


79 

m 

79 

.000 

c 




1 


n 
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So 

m 

SO 

*000 



IF (CnttGT *2(3. TGTID) RETURN 


ai 

m 

81 

*000 



NXTHDG » TGTHDG(CnTTGT) 

82 

m 

32 

*000 



IF OLOHOG *GT. NXTHDG) NXTMDG » NXTHDG + 

TW8=I 

S3 

9* 

83 

*000 



oldhog • nxthdg 


84 

m 

84 

*000 



If (NEWHOq *gT* NxTHDG) g9 T6 30 


85 

m 

85 

*000 

c 




86 

M 

36 

*000 


CHECK IF MOVE NECESSARY* •* 


87 

m 

87 

*000 

c 




88 

■i 

88 

*000 



If (lentrv *ne* 1) g9 re 35 


89 


89 

*000 



RETURN 


90 

m 

90 

*000 

c 




91 

m 

91 

*000 


REMOVE THE TARGET*** 


98 

w 

92 

.000 

c 




93 

m 

93 

*000 


35 

IDCCW a TGTPTR(nTQTID) 


94 


94 

*000 



IDCW = T3TPTR(2/ TQTID) 


95 


95 

.000 



TGTPTR(2, IDCCW) = lOCW 


96 

m 

96 

*000 



TGTpTRd.IDCW ) * IDCCw 


97 

m ' 

97 

*000 

c 




9S 

m 

98 

*000 


INSERT THE TARGET.** 


99 

m 

99 

.000 

c 




100 

m 

100 

*000 



IDCCW s TQTPTr( nCNTTGT) 


101 

nt 

101 

*000 



TGTPTR(2^ lOCC'rt) s TGTID 


102 

n 

102 

*000 



TGTPTR(1#TGTID) » IDCCW 


103 

m 

103 

.000 



TGTPTf’(2/TGTlO) s CNTTGT 


lOii 

«* 

104 

.000 



TGTPTR( 1#cNTTGT) a TGTID 


105 

m 

105 

*000 



RETURN 


106 

m> 

106 

*000 

c 




107 

m 

107 

.000 


traverse counterclockwise* • • 


108 


103 

*000 

c 




109 

m 

109 

*000 


*^0 

CNTTGT a TGTID 


no 

m 

no 

.000 



LENTRV a 0 


111 

m 

111 

*000 

c 




112 

n 

112 

*000 


45 

LENTRV = LENTRv + 1 


113 

n 

113 

• 000 



CNTTGT = TGTPTR{ 1/C.NTTGT) 


114 

m 

114 

.000 

c 



115 

m 

115 

.000 



If (CNTTGT *EQ* TqTID) RETyRN 


116 


116 

.000 

c 




117 


117 

*000 



NXTHDG a TGTHDG(CNTTGT) 


118 


113 

*000 

c 




119 


119 

*000 



IF (OLOHOG *LT* NxTHDG) NxTHDG = NxTHOG - 

TW6PI 

I2n 


120 

•000 



SLONDG a nxthdg 



Ill 



121 

m 

121 

•000 

IF (NEWHOg •LT. NxTHOG) g9 T8 45 

—1 

12? 

m 

122 

• 000 

c 

j 

12.1 

w 

123 

• 000 

c*** CHECK IF M 0 VF NECESSARY... 

i 

124 

m 

124 

.000 

C 


125 


125 

• 000 

If (iENTRV .NE. 1) gS T8 50 
rfturn 

n 

i 

126 


126 

• 000 

. ] 

127 

• 

127 

• 000 

c 


12.^5 


12S 

• 00 c 

c*** RFMOVE the target... 

n 

139 

m 

139 

.000 

c 

1 

j 

13Q 

m 

130 

• 000 

50 loccw * TGTPTRC li TQTID) 


131 

m 

131 

• 000 

lOCW s TGTPTR(2>TGTI0) 


13? 

«■ 

132 

*000 

TGTPTRf2i IDCCW) « IDCW 

1 

133 

m 

133 

.000 

TGtPTR( 1. lOCW ) 3 lOCCW 

. 1 

134 

135 

m 

134 

135 

• 000 
• 000 

8*** TNSFRT the TaRqET... 

n 

1 

136 

m 

136 

• 000 

c 

1 

. ) 

137 

m 

137 

• 000 

lOCW » TGTPTR(2/CNTTGT) 


13a 

m 

138 

,000 

TGTPTR( 1, lOCW) a tGTIO 


139 

m 

139 

• 000 

TgTpTRf 1 ,TqTI 0) » CNTTnT 

I 

140 

m 

140 

.000 

TGTPTR(2.TGTID) >? IDCW 
TGTPTR(2,CNTTGT) = tgtio 

. 1 

141 

m 

141 

,000 


14? 

m 

142 

.000 

return 

n 

143 

m 

143 

.000 

c 

. 1 

144 

m 

144 

.000 



145 


145 

.000 

END 9f subroutine »ScIUpD'... 


146 

- 

146 

• 000 


l47 

m 

I 47 

.000 

c 

. 

148 

m 

148 

• 000 

ENq 



1 

'n 
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SUBROUTINE HDGDST 


DESCRIPTION: 

COMPUTES HEADING AND DISTANCE OF GIVEN TARGET USING 

THE DISTANCE OF THE HELICOPTER FROM THE PRIMARY TARGET AND 

THE DISTANCE OF THE TARGET FROM THE PRIMARY TARGET 

COMMON DESCRIPTION: 

SEE SUBROUTINE CREATE. 

VARIABLE DESCRIPTION: 

DELX - DIFFERENCE IN X LOCATION 

DELY - DIFFERENCE IN Y LOCATION 

FT2NM - FEET TO NAUTICAL MILE CONVERSION 

NEWHDG - NEW HEADING OF TARGET 

TGTID - ID OF TARGET BEING UPDATED 

SUBROUTINE CALLED: 

NONE. 
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NAMf = 

1 • 

a 

2 - 

4 w 

5 * 

ft ^ 

7 - 

8 - 
9 - 

10 - 
11 - 
12 - 

1 3 " 

14 - 

15 - 

16 * 

17 - 

18 • 

19 - 

20 - 
21 • 
22 - 

23 - 

24 - 

25 , 
2ft * 
27 ■» 
23 - 
29 - 
3Q ^ 

31 - 

32 ^ 

33 • 

34 - 

35 - 

36 - 

37 - 
33 - 


S:H0G0ST.HEL9An 


1 

2 

3 

4 

5 

6 

7 

S 

9 

10 

11 

12 

1 3 

14 

15 

16 
17 
IS 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


000 

00c 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 


c*# ***»***♦■*■*»*•** **#*4* *■**■**♦♦ *4*4*4 •»*»■* ***»*##* 4 **lt#*» 

c**************** 4 *** 4 *** 4 ** 44 * 4 ** 44444 * 44 ** 4 ** 4 *>*44*444*4**** 

c***** TITi_E..» S;HqGq 3 T> HEi ****** 

c***** last modified; 20 jaN/ "as: vdhe.'j 4 ***** . 

c*** THIS subroutine computes the heading and distance of 

C**4 GIVEN targets USING THE X- AND Y- COORDINATES OF THE — , 

c*** Target an^ helIcSpter. . f no EllIpticAl csppEcTions I 
c*** are used*** 
c*** 

c*************************************** ************** ********* I 

c^** SuBr 8 uTI-'^ES called 5 NONE**. ! 

C*******#*******»*******#4**4**4******4********4**4************ 

c - 

SUBROUTINE HDSDST 

c . 

c************************************************************** 

c***4 * common blocks*** **♦**#! 

C**************************************************************' I 

c 

COMMON 

* /XFLOAT/ A(500 )/IFIXED/I a ( 250) I 

* /HRCeM /RH( 50), IH(50),XC( 2b) 

COmm®'^ — ' 

* /BINTRE/9INLST(5i 150)/ HDBLST, NXTN0D i 

* /TARQET/TGTLQC(2i lOO)# TgTDST(IOO)/ TQTHDQ(lOO)f ‘ 

* TGTPTR(2i lOOji TGTLOO(IOO)^ TGTTLT(IOO), 

* tGttyP(Ioo), n 

* NBRTGTi NysiTGT/ TGTID ; 

* /STATUS/LSTDIRi NxTAlT/ NXTUPD/ TGTEnD, TGTSTRi LaZM 

C -1 

»*****•»’»♦*'*****♦* »**♦*♦** »■>!>#*«'*'» »*»***** 
C***#* EqU I VALENCES • #♦♦**♦ * 


Equivalence 
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— 3s - 

39 

.000 


( A{103)j XPR)i ( A(104)i YPR)i 1 A(105)i HPr), 

40 • 

40 

*000 


{ A(353)i D2R)i ( A{460)i PI)/ ( A(461)/ TWBPT) 

41 , 

41 

.000 

c 


42 - 

42 

.000 


^ 43 

43 

«00C 

c ♦•»*»* * 

DECLARAT I dNS • • • **##<»■# 

44 - 

44 

.000 


45 - 

"♦5 

.000 

c 


. — 46 •> 

46 

.000 


REAL NEwHDG 

47 • 

47 

.000 

c 


4s 

48 

• 000 


IMTEGER TGTEND# TGTIO# TuTPTRi TGISTR, rGTTYP 

49 - 

49 

.000 

c 


50 - 

50 

.000 


5t - 

51 

*000 


DATA STATEMENTS... ****** 

52 - 

52 

•000 


— 53 - 

53 

.000 

c 


54 ^ 

54 

• 000 


qATA FT2NM/1 .6447363E-4/i XPp/- 12 1 60 . 0/i YPr/O.O/ 

55 - 

55 

.000 

c 

_ 56 - 

56 

.000 


57 - 

57 

• 000 


5fl - 

58 

.000 

c 


59 - 

59 

.000 

c* **** 

ceMPuTE The distance... 

- 60 - 

60 

.000 

c 


61 « 

61 

.000 


OElx = TGTl9C{1/TGTId) - XPr 

62 " 

62 

.000 


DELY = TGTL0C(2^TGTID) - YPR 

_ 63 •* 

63 

•000 

c 

6 

64 

.000 


TGTDST(TGTID) = SQRT(0ELX»»2 + DELY**2 + HPR**2 ) *FT2NM 

65 m 

65 

,000 


TGTL3G(tGtID) » L0G1O(tGT'DST(tGTID) ) 

66 •• 

66 

*000 


ALTNM » -HpR*FT2NM 

57" 

67 

.000 


TGTTLT(TGTID) « ASIN( ALTNM/TGTDST(TGTIDn 

63" 

68 

.000 

c 


69 - 

69 

.000 


IF (DELx .Eq. 0.0) RETURN 

- 7o - 

70 

.000 

c 


! 71 - 

71 

.000 


COMPUTE The heading,, , 

^ 72 " 

72 

• 000 

c 


73 " 

73 

.000 


NEWHOQ » ARCrAN(DELY/DELX) 

~ 74 - 

74 

• 000 


IF (NEWHOG .lT, 0*0) NEWHDG = NEWhDG + Tw8PI 

75 " 

75 

,000 

c 

76 •’ 

76 

.000 


T f5 ThD G ( Tq T I D ) s NEwHDg 

— 77 ^ 

77 

•000 

c 

» 

73 " 

73 

*000 


return 

79 - 

79 

.000 

c 
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80 

m 

80 

*000 


. 81 


81 

• 000 


82 

m 

82 

*000 

c 




83 

m 

83 

.000 


RFAL FUNCTieN ARCTAN ( DX# DY ) 

34 

m 

34 

,000 

c 




35 

m 

85 

*000 


DATA HALFpI/1.5707953/i pI/3. 14lb9P7// T»0pI/6. 2331353/ ~ 

86 

w 

86 

*000 

c 




37 

m 

87 

*000 


IF 

( aBS(OX) 

•Gr. ABS(DY)) G 0 T8 iQ 

88 

W 

38 

.000 



IF (DY 

.£q. 0.0) ArcTAm 3 0«0; RETUR^J -] 

89 

m 

89 

.000 

c 



i 

90 

m 

90 

.000 




ApCTAN = ATAN(qX/dy) ‘ 

91 

m 

91 

.000 


IF 

,DY 

•LT. 0.0) ARCTAN _ ARCfAN + PI 

9? 

«• 

92 

•000 


IF 

{ arctan 

• GT. PI) ARCTAN = ARCTAN - T.vOPi j 

93 

m 

93 

,000 


RfjTURN 

. 1 

94 

m 

94 

.000 

c 




95 

m 

95 

.000 

10 



A-RCTAN = HALFPI - ATAN(OY/Dx) -; 

96 

w 

96 

.000 


IF 

(OX .LT* 

Q.QJ ARCTAN » ARCTAN - PI j 

97 

m 

97 

.000 


RrTURN 


98 

m 

93 

•000 

c 



•— 1 

99 

m 

99 

.000 


100 

■» 

100 

.000 

c* 


end 9F SUdRsUTINE ’HDGDST’##* *#***#, 1 

101 


101 

.000 


102 

*• 

102 

.000 

c 



”1 

103 

m 

103 

,000 


E\'D 


1 


"1 

i 


I 


“1 

i 
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SUBROUTINE INIROR 
DESCRIPTION: 

INIROR INITIALIZES THE TARGET DATA NECESSARY FOR THE 
RADAR SIMULATION. THE HELICOPTER LOCATION MUST BE 
INITIALIZED IN COMMON XFLOAT PRIOR TO CALLING THIS 
ROUTINE. ALL THE TARGETS ARE DATA INITIALIZED IN COMMON 
/TARGET/. IN ORDER TO CHOOSE THE BLOCK OF TARGETS. THE 
OPERATOR MUST FIRST SET NURTGT TO THE FIRST TARGET IN 
THE BLOCK. THEN A CALL TO INIROR WILL INITIALIZE THE 
TARGET DATA FOR THAT BLOCK OF TARGETS. 


COMMON DESCRIPTION: 

SEE SUBROUTINE CREATE. 

VARIABLE DESCRIPTION: 

NONE. 

SUBROUTINE CALLED: ‘ 

INSTGT - INSTALL TARGETS INTO TARGET LIST. 



name s 


I - 
? - 
3 - 

k - 

5 - 

6 - 
7 - 

a - 

9 - 

10 - 
11 - 
12 - 
13 - 
U - 

15 " 

16 - 
1 ? - 
13 - 

19 - 

20 - 
21 , 
22 - 

23 - 

24 - 

25 - 

26 - 
27 - 
23 - 

29 - 

30 - 

31 - 

32 - 

33 - 
3 ^ . 

35 *■ 

36 « 

37 - 

38 


SJ INIROR .HELRAD 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
23 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

ooc 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 


^•»#******>*#*****# »**♦#*****»******♦»»■»*»**»**» »»»»*****»»**4** 

C*-»****-*-*******»**»***'»**»»***»»******»****»****** *##»♦**•***#*•* — : 

c***** title.** S:Inir[dr/ helrao **#*♦«> : 

LAST MQOIFIED: 20 JANj ’s2; VDbEN ^,*^#**j 

^M^*#*****#^MHHM|.#*#*l^**♦*^^#*^MH^**^M^*^^•)MM^♦4******R^**^^*#**»*^M^#^f4^^ 

c*** 

C*** THIS SuBRauriNE initializes The target data necessary 
c*** FOr the HElIcQ^TEr rAqAr SIMULATION... AlL TArGETS AdE 
c*** inserted vIA a call T8 'INSTGT'..* the parameters are 

C*»* THE X- AND Y- C88RDINATES 6F THE TARGET (iN FEET)... 

C*** THE HELlCaPTER SIMULATIBN MUST 3E INITIALIZED IN 

C*«# IxFlBAT* C3MM9N PRI3R T8 CALLING THIS SUBROUTINE... H 

C**** #♦♦**»♦*#•»***** ****#■»* ♦♦♦*9** #**♦*-* ♦#**♦#'♦* 

C*** SUBROUTINES CALLED: 1 


in'STqT - install tarqet data iNra the tarqet list... 


1 


c**** 

c * * * 1 

C********#***-***** )*♦***»* *•»*♦****♦**•»***»******♦*♦♦*•**» »**♦*♦*« I 

c 

SUBROUTINE INIRDR 
C 

C#**** *****♦♦♦**■****»■* *■»*»***♦*»*•**»■**»***♦***** **♦*•* **!.** **♦*■* 
C#**** COMMON SLOCKS*.. **#*♦*-. 

C*»***** *»•»■*********»♦♦»***»*♦*»•********»****■•**»**»**» ****»**•« ; 

C3MM8N 

/XFL8AT/ A(500) / 1 F I XEO/ I A ( 25o ) i 

/HRCaM / RH(50), iH(50), KC(25) 

C^MMaN 

/BINTRE/aiNLST(5/150)i HDBLST, NXTNao 

/TARgET/TgTl8C{2/100). TgTDST(IOO)/ TGTHDG(IOO). j 

TGTPTR(2/ 100)# TGTlOG(IOO)# TGTTlT(IOO)# 

TGTTYP{100), _ 

NBRTGT# NUMTGT# TGTID ! 

/STATUS/lSTdIrj NxTAlT# nxTuPdJ' TQTENo# TGTSTr# laZP i 


000 c 
000 
000 

1 

000 

* 

000 

000 

1 

# 

000 

* 

000 

* 

000 

«• 

000 

* 

000 

# 


n 

1 
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I 


3q ^ 

40 - 

41 - 

42 - 

43 - 

44 " 

45 • 

46 * 

47 - 

48 - 

49 •• 

5 0 ^ 

51 - 

52 ^ 
5:^ » 

54 • 

55 • 

56 !• 

57 - 
5S ' 

59 ^ 

60 - 
6 1 *• 

6 ? • 

63 " 

64 ” 

65 

6b - 

67 •> 

6i^ 

69 ^ 

70 - 

71 - 
7? - 

73 - 

74 - 

75 -r 

76 - 

77 - 

78 - 
73 - 


39 

40 

41 

42 

43 

44 

45 

46 

47 

43 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 
63 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 


/ANTE^IA/ANJTOIMi RAOdi-lEj ANGAiN, BRKP.nTj RNG3R:</ AMTL8G 


000 * 

000 C 

000 C* *♦**•» »* 4 ******** ******* *******♦*#)►*»**»****♦»•»»•♦•»** *•**♦»*•**•*■» 
000 c***** OECLARATISNS . • . ****** 

ooo c******** ********************************************** ******** 

000 c 

000 INTEGER TGtEND*. T 3 TID/ TGTPTR/ FGTSTRi fGTTYP 

000 C 

000 c ****************************************************** ******** 

000 C***** DATA STATEMENTSt » « *«**»» 

000 c******** ********************************************* ********** 

000 c 

000 DATA ANGAIN/ 28 » 0 /i ANTD I M/ 1 2 . 0 /> TGTPTR/ 200 * 0 /i 

000 * X'^GS/ 500 . 0 /, XRL 9 C/ 1500 ,/i 

000 * Yf^GS/-100«/i yRLdc^O'OOO/ 

000 C 

000 data NRRTGT/ 3 /> NUMTGT/ 7 /i TGTTyP( 3 )/ 1 / 

000 C 

000 C *'**********♦*******»***♦**********♦**** **♦****♦*^^******♦**^^*** 

ooc c***** target data#** ****** 

000 c************************************************************** 

000 c 

000 C*** ILS RADAR targets*** SET »NBRT 3 T» = 1 ; f.sjUMTGTt = 2 ... 
000 C 


*000 


DATA 

TGTL9C( 1^ 

1)/ 

500*0/:» 

TgTL8C(2, 

1)/ 

-100*0// 

•ooo 

,000 

c 


TGTl9C(1> 

2)/ 

1500*0// 

T3TLeC(2i 

2)/ 

0.0/ 

• ooo 
•ooo 

c* ** 
c 

6IL RIg 

targets.* . 

S£ 

:t ‘nbrtgt 

' » 3; »nU,v(TgTi 

= 7.** 

• ooo 


DATA 

TGTl9C{ 1* 

3)/ 

0*0// 

TGTL9C(2* 

3)/ 

0*0// 

• ooo 


* 

TGTL0C(li 

4)/ 

1986*0// 

T3TL8C(2, 

4)/ 

-304*0// 

• ooo 



TGTL9C(d 

O)/ 

3972*0// 

TGTL0C(2/ 

5)/ 

2723*0// 

,000 


* 

TGTLGCd, 

6)/ 

3986*0// 

TGTL9C{2* 

6)/ 

561.0// 

• ooo 


■» 

TGTLQCd# 

7)/ 

3669 * q// 

TGTLQC(2i 

7)/ 

-1430*0// 

• ooo 


* 

TGTLQCd* 

8)/ 

992*0// 

TGTL3C(2/ 

8)/ 

9)/ 

-3491*0// 

,000 
• ooo 

c 

♦ 

TGTL9Cdi 

9)/* 

80g4,0// 

TGTL8C(2, 

-5610,0/ 


000 C*** SANTA BARBARA CHANNEL. (9 TARGETS) ••* 

ooo c*** 

000 c 

ooo DATA 


Set 'NBRtGT’ * lO; 'NUMrQl' » 9 *.. 
TGTLSCdi 10)/-2l400*0/i TGTLSC(2* 10 ) /- 1000* O/i 


119 



Sc 

«■ 

30 

•000 


* 



TGTLaCf u 

11)/ 

^13500. 0/i 

TGTLSCf?/ 

11)/ ■9oo.n// 


81 


31 

• 000 





TGTlOCI li 

12)/ 

■15800. 0/^ 

rGTL8C( 2/ 

12)/ ■200.0// 


S3 

m 

82 

•000 


* 



TGTuaCf 

13)/ 

■13100.0// 

TgTL8C( ?/ 

13)/ -160.0// 


S3 

« 

S3 

*000 


♦ 



TGTl9C( 1» 

l4)/ 

0*0/* 

TGTL8C(?/ 

1^57 0*o7* 


34 


84 

,000 


* 



TGTL3C( 1, 

15 ) / 

2380,0/, 

TGTL8C(2, 

'15)/ ,g50, 07, 

-0 

35 


35 

•000 


* 



TclLacC 1/ 

16)/ 

4990*07/ 

TqTL8C(2/ 

16)7-1 3nO. o7/ 


86 

m 

86 

.000 


4k 



TGTL3C( 1# 

17,/ 

5480.0// 

TGTL8C(2/ 

17)/-2400.0// 

( 

«7 


8^ 

,000 


4k 



TGTLaCI 1, 

Igi/ 

12120,0/, 

TGTL8C(2/ 

lg)/-3200,07 


85 

a* 

38 

• 000 

c 









; 

33 

W 

39 

• 000 



1 

nautical mile arc... 

SET iNBRTGT' = 19; 

i.NUMTGTi 3 3.. 

1 

. I 

gO 

if 

gO 

,000 

c 










91 

m 

91 

•000 


D 

ATA 

TGT|_8c ( 

19)/ 

5364.9// 

TGTl8C(2/ 

19)7 -2852*67/ 


92 

m 

92 

• 000 


* 



TGTuaC( J / 

20)/ 

5709.7// 

T3TL8C(2/ 

20)/ -2078.1// 

. j 

93 

m 

93 

•000 


4k 



TGTl3C( 1/ 

21)/ 

5943*3// 

TgTL9C<2/ 

21)/ -1263.3// 


94 

m 

94 

•000 


4k 



TGTl8C{ 1/ 

22)/ 
23) / 

6q61.3// 

TQTL8C(2, 

22) 7 -A23.8// 

23) / ^23.,/, 

24) / 1263.3// 

n 

1 

q 5 

if 

o5 

,000 


4k 



TGTLaCd, 

6061,3// 

TGTL3C(2, 

96 

m 

96 

•000 


4k 



TfjLaCf li 

24)/ 

5943.3// 

T3TL8C(2/ 

- i 

97 

m 

97 

•000 


4k 



TGTL3C( 

25)/ 

57Q9 * 7_// 

TGTL8C(2, 

25) 7 2078.17/ 

26) 7 2h52.67 


93 

m 

93 

,000 


4k 



TGTL3C( 1, 

26)/ 

5364, g/. 

TaTL0C(2, 

“] 

99 

m 

99 

• 000 

c 








1 

100 

if 

100 

.000 

c*** 


2 

nautical mile arc... 

SET tNBRTGT* = ? 7 i 

iNUiMTGTi s 3.. 

• 

lot 

«» 

101 

.000 

c 










10? 

m 

102 

• 000 


D 

ATA 

TGTlQc ( 1 * 

27)/10729.8/j 

^GTl8c(2/ 

27)/ -5705.1// 

1 

103 

m 

103 

• 000 


4k 



TGTl6C( 1/ 

23)/114l9.3// 

TGTL8C(2/ 

28)/ -4156.3// 

I 

1 

10^ 

m 

10^ 

•000 


4k 



TGTl3C( 1/ 

29)/lia86.6/^ 

TGTL8C(2/ 

29)/ -2526.67/ 


105 

m 

105 

,000 


» 



TGTl8C( 1/ 

30)/12l22,6// 

TQTL8C(2/ 

30 ) 7 - 347 .77/ 

i 

106 

m 

106 

•000 


♦ 



TcTuaCf 1/ 

3D/12122.6// 

TqTL9C(2/ 

31)/ 847.7// 

1 

- ] 

107 

m 

107 

• 000 


* 



TGTLOCdi 

32>/11386.6// 

TGTL0C(2, 

32) 7 2526.6// 

33) 7 4156,37/ 


10^ 

m 

108 

,000 


4k 



TGTLacd, 

33)/ll41g,3/, 

TGtL0C(2, 

n 

103 

m 

109 

• 000 


4k 



TgTl 3C( li 

34)/10729.?/i 

TgTL8C(2/ 

34)7 5705.17 

1 

110 

m 

110 

,000 

c 









. J 

111 

m 

111 

.000 

c*»* 


3 

nautical mile arc... 

SET INBRTGT* » 35; 

»NUMTGT» » 8.. 

• 

112 

m 

112 

.000 

c 









[ 

113 

if 

113 

.000 


DATA 

TGTL8C( li 

35)/16094.6/i 

TGTL8C<?/ 

35)7 -8557.77/ 

\ 

114 


114 

• 000 


4k 



TGTl8C(1^ 

36)/l7l29.0/^ 

TGTL8C(2/ 

36)7 -6234,47/ 


115 

•* 

115 

• 000 


4k 



TGTl8C(1^ 

37)/17830.0// 

TGTL0C{?/ 

37)/ -3739.9// 

— 

116 


116 

•000 


4k 



TGTl8C( 1/ 

38)/13l83.9// 

TGTLaC(2/ 

38)/ -1271.57/ 

i 

1 

117 

m 

117 

.000 


4k 



TGTL8C(1. 

39)/13lS3.9/i 

T3TL8C(2/ 

39)7 1271,57/ 

I 

ils 

m 

11a 

.000 


4k 



TGTl 3C{ d 

40)/l7830,0/i 

TGTL8C(2/ 

40)7 3789.97/ 


119 

m 

119 

.000 


4k 



TGTl8(^ ( 1 * 

41 )/17129.0// 

TgTL8c(2/ 

41)7 6234.4// 

"1 

120 

m 

120 

.000 


4k 



TGTL3C( 1/ 

42)/16094.6/^ 

T3TL8C(2/ 

42)/ 3557.7/ 

1 


'n 
, 1 


120 



t* MAUTICAL mile arc**. 


SET «NBR TQT' 


43; •NU'^TGT' a a**. 



' 121 

m 

121 

*000 

c 


1 122 

m 

122 

*000 

c *** 


■'123 


123 

.000 

c 


r -124 

w 

124 

*000 



i 125 

m 

125 

.000 


♦ 

■ 1?6 

m 

126 

.000 


- M - 

127 

m 

127 

*000 


* 

r 123 

■■ 

128 

*000 


* 

i 129 

•• 

129 

*000 


* 

130 

m 

130 

*000 



-131 


131 

,000 



132 

■t 

132 

• 000 

r 


133 

■1 

133 

•000 

c »« • 


^34 

** 

134 

•000 

c 


: 135 


135 

•000 



136 

wm 

136 

,000 


* 

137 

m 

137 

• 000 


* 

;“138 

m 

138 

*000 



i I 3 g 

m 

139 

*000 



l 40 

m 

140 

.000 



-141 

m 

141 

♦ 000 


* 

. 142 

m 

142 

*000 



143 

m 

143 

•000 

c 


144 

m 

144 

• 000 

c *» ♦ 


"145 

mk 

145 

.000 

c 


146 

m 

146 

*000 



147 


147 

*000 


* 


m 

143 

.000 


* 

149 

m 

149 

*000 



i 5 o 

m 

150 

*000 



151 

W 

151 

.000 


* 

“^52 

m 

152 

*000 


* 

153 


153 

• 000 



154 


154 

*000 

c 


—155 

m 

155 

•000 

c*** 


; 156 

m 

156 

.000 

c*-** 


157 

n 

157 

*000 

c 


153 

m 

158 

•000 



. 159 


159 

.000 



I60 

W 

160 

• 000 


* 

161 

m 

161 

*000 




DATA T3Tl9c(1» 43)/2l459*5/i T3Ti_8c(2i 43') /-11410. ?// 
TGTl 8 C( 1 / 44)/2283S*7/^ TGTL9C(2i 44)/ -8312. A/i 

TQTLaCdi 45)/23773,3/^ ri3TL3C(2^ 45)/ -5053.?// 

TGTL^cd/ 46)/24245.2/i TGTL8c(2i 46)/ -1695*4// 

TGTl 9C(1/ 47)/24245.c// TGTL9C(2, 47)/ 1695*4// 

TQTl 3C(1/ 4S)/23773.3// T3TL9C(2i 48)/ 5053*?// 

T 3 Tl 9 C( 1 / 49)/2283a*7// T 3 TL 6 C( 2 * 49 )/ 8312 * 6 // 

TGTLdCd/ 50)/21459*5// TQTL0C(2. 50)/ 11410,2/ 

5 NAUTICAL MiLE ARC*.. SET iNb^ TqT i s 51; ’NUMTqT' » 

DATA TGTLSCd/ 5 1 ) /26S24 * 4// TGTL0C(2/ 51 ) /- 1 4262 * 8 // 
TGTLQCd, 52)/2S548,3/, TGTL 0 C( 2 . 52 ) /-lOSgO , 7 /, 

ToTLaCd/ 53)/29716.6// TqTL 8 C( 2 / 53)/ -6316*5// 

TGTL0Cd/ 54)/30306*5// TGTL9C(2/ 54)/ -2119.?// 

TGTLdCd/ 55)/30306.5// TGTL9C(2* 55)/ 2119,?// 

TQTl9C(1/ 56)/29716*6// TqTL9C(2/ 56)/ 6316*5// 

TGTLOCdi 57)/23548.3// TQTl9C(2/ 57)/ 10390*7// 

TGTLecd^ 53)/26824.4// TGTL9C{2* 58)/ 14262*8/ 

10 NAUTICAL l-llLE ARC*** SET »N3R TGT i s 59 ; iNUMTGTt * 

DATA TGTl9C(1/ 59)/53649*0// TgTL8C(2* 59) /-28526.Q// 
TGTl9C(1/ 60)/57097.0// T3TL0C(2* 60 j /-20781 , 0 // 

TGTLQCd, 61)/5g433,0/> TGtL9C(2, 61 )/- 12633, 0 /^ 

TQTL9Cd/ 62)/60613*0// TqTL9C(2/ 62)/ -4238*0// 

TGTLeCd/ 63) /60613.0// TGTLeCI?/ 63)/ 4238*0// 

TGTL9C(d 64)/5g433.C/, TGTL9C(2. 64)/ 12633,0// 

TGTL9C(1/ 65)/57o97.0// TgTL9C(2. 65)/ 20781*0// 

TGTL9Cd/ 66)/53649*0// TQTLeC( 2 * 66 )/ 23526 * 0 / 

SANTA BARBARA CHANNEL (11 TARGETS),., 

set 'NBR tGT' = 67; *NUMtGt' » 

DATA TGTLSCd/ 67 ) /-2l 400 * 0 // TGTL 9 C( 2 / 67) /-looo*o/.» 

TGTL9C(1/ 63)/-l3500,0// TGTL 8 C( 2 / 63)/ ,900,0// 

TgTL9C(1/ 69)/-15800.0// TgTL0C(2/ 69)/ -200*0// 

TGTLaCd/ 70)/-13100.0// TGTLBC(?/ 70 ) /' -I6O.0// 


8 * * ♦ 


8 • * • 


11 * • 


121 



162 

i* 

162 

• 000 

163 

m» 

163 

•000 

164 

m 

164 

.000 

165 

m 

165 

• 000 

166 

m 

156 

•oon 

167 

«■ 

167 

.000 

163 


158 

• 000 

169 

m 

169 

• 000 

170 

m 

170 

• 000 

I7I 

tm 

I7I 

.000 

172 

m 

172 

.000 

173 

m 

173 

• 000 

174 

m 

174 

• 000 

175 

w 

I75 

.000 

176 

w 

176 

• 000 

177 

■■ 

177 

• 000 

173 

w 

178 

• 000 

179 

m 

179 

• 000 

I80 

W 

180 

• 000 

I8i 

m 

181 

• 000 

lfl2 

w 

132 

,000 

183 

m 

183 

• 000 

184 

m 

184 

• 000 

185 

m 

185 

*000 

186 

m 

136 

.000 

187 

w 

187 

• 000 

138 


188 

.000 

189 

m 

189 

.000 

1^0 

w 

190 

•000 

191 

m 

191 

• 000 

192 

m 

192 

.000 

193 

m 

193 

• 000 

194 

m 

194 

• 000 
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LSTDIR = 
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LAZm a 
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ANTL 9 G = UeGlO(ANTOlM) 

BRKpNj » 10,0*»((,2.4 + ANqAIN)/16.86) 
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RETURN 
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END 


END 0 F SUBR 8 UTINE MNIRDRi... 


*■»#*** 
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SUBROUTINES INSTGT 


DESCRIPTION: 

INSTGT INSERTS INTO THE TARGET COMMON AREA THE COMPUTED 
BEARING AND DISTANCE FROM THE HELICOPTER AND MAINTAINS 
THE TARGET POINTERS SO THAT THE TARGETS ARE IN A CIRCU- 
LAR LIST IN ORDER OF BEARING. ALSO DEFINED IS THE 
STARTING VALUES OF VARIOUS POINTERS. INTO THE TARGET 
LIST. 


COMMON DESCRIPTION: 

SEE SUBROUTINE CREATE. 


VARIABLE DESCRIPTION: 

FSTVIS - FIRST VISIBLE TARGET IN TARGET LIST 
HOGOIF - HEADING DIFFERENCE BETWEEN TARGET AND HELICOPTER 
HDGEND - HEADING DIFFERENCE OF LAST VISIBLE TARGET 
HDGSTR - HEADING DIFFERENCE OF FIRST VISIBLE TARGET 
HEAD - HEAD OF TARGET LIST (SMALLEST HEADING) 

IDCCW - COUNTERCLOCKWISE TARGET 1,0. 

IDCW - CLOCKWISE TARGET I.O, 

LSTVIS - TARGET I.O. OF LAST VISIBLE TARGET 
NTGHDG - HEADING OF TARGET BEING INSERTED 
TAIL - TAIL OF TARGET LIST (LARGEST NUMBER) 

TGTID - I.D. OF TARGET BEING INSERTED 
XLOC - X LOCATION OF TARGET 
YLOC - Y LOCATION OF TARGET 


SUBROUTINE CALLED: 

HDGDST - COMPUTES HEADING AND DISTANCE OF TARGET 
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c***** titue».» S; 

last mSoIfIED! 20 JA 


c*** This subRQutine inserts iNie the target cqmhS^ area the 

C*** X- and Y- ceSf^DlNATES BF ThE TArGETS; CQmPuTES T^E BEAr- 
C*** ING AND distance FR9M THE HELICSPTeR AND MAINTAINS THe 
c*** target pointers se that ThE targets are in a circular 

C*** LIST IN ORDER OF THE BEARING yALUES... 
c** * 


S.lNSTGTi HELRAD 
20 JAN, ' 32 ; VDbEN 






c** * 
c*-* * 
c*#* 


SU^ReUTlNES CALLED: 

HOQDST - CSMPUTES THE HEADING AND DISTANCE QP EACH 
T ArgET • • • 


Subroutine instgt 


C****» CSMMON 3LBCKS».t **#♦** 

C********)MHHHH»***** ******•»••*♦*•*»*#♦*♦**»»**»**#* ****»*■*»#♦****♦ 

c 

common 

* /XFLOAT/ A(500)/IFIXED/IA(250) 

* /HRCqM /RH( 50)^ IH(50),KC( 2b) 

COmmBN 

* /BINTRE/BlNLST(5il50), HDBLST, NXTNSQ 

* /TARgET/TgTl 9C(2^ lOO)/ TgTDST(iOO)^ TGTHDG(iOO)> 

* TGTPTr(2/ 100)# TGTl9^(100)# tGtTlT(IOO)# 

* TGTTYP(IOO), 

* N3RTGT, NUMTGT# TGTIO 

* /STATUS/LSTOIR# NxTAlT, NXTUPD# TGTEnD# TGTSTR, LAZM 


****** 
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C * * ♦ 
C**i 
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« 
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DECLARATI6NS. . . 


REAL NTgHDG 
integer 


FSTVIS> HEADi TAlLi TGTEND, 
TGTID/ TGTPTRi TGTSTR^ TGTTYP 


C* * ■)*■# * 
C*#-*##****-! 

c 

data 

* 

c 


data statements*.. 


•»•*#•*.** 


Pl/3.lA15927/i TWaPI/S. 2831353// mAxTGT/100/ 


c 

c*** get new target id... 


c 

IE (NUMTGT .GT. (MA/TGT-NbPTGT) ) NUMTGt = MAXTGT-NbRTGT 

c 

IE MNBRTrjT ,LE. 0) .6R. { NbRTqT+NUMTqT- 1 .gT. MAxTqT)) 

* RETURN 
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INITIALIZE THE TArGET PQINTEr A?^rAY TB ZEr0**. 
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TqTpTR ( I / J ) 3 0 
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random number and produces a new integer And real RANn0M_^ 

N’UMBEr ... j 


1 

n 
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peal v,|rLp 
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iNT-fitR INTx> INTY 
C 

INTY = I'^TX*65539 
iNJx » IMTY 
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c 
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IE (INTY)/ 1/1/2 


1 
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INTY = INTY + 2U7483647 + 1 
WFLP = INTY 

WFLP a WELP*0.4656613E-9 
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return 

c 

C***** end 0F subroutine 'RANDUMI.,, ****♦» 

C 


END 
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